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SECTION  I 
INTRODUCTION 


1.  OBJECTIVES 

In  designing  the  Air  Force  Refuse-Collection  Scheduling  Program  (RCSP), 
the  fundamental  objective  was  to  reduce  collection  costs.  The  most  signifi- 
cant cost  reduction  Is  effected  by  a reduction  In  the  number  of  collection 
trips  used  to  service  a given  region.  If  a collection  crew  can  be  dropped 
from  the  fleet,  the  cost  of  manpower  will  be  cut.  In  addition,  fuel  and 
maintenance  costs  will  be  lessened  If  the  total  mileage  traveled  by  the  col- 
lection fleet  can  be  reduced.  The  first  objective,  then,  Is  to  generate  a 
collection  schedule  that  calls  for  the  theoretical  minimum  number  of  trips. 
This  objective  Is  accomplished  In  program  PHASE2  of  the  RCSP. 

A secondary  objective,  good  spatial  clustering  of  all  streets  serviced 
by  a vehicle  during  one  trip,  Is  also  achieved  by  PHASE2,  except  possibly 
for  the  last  trip.  In  addition  PHASE2  plots  maps  that  show  the  section  (trip) 
to  which  each  street  segment  Is  assigned. 


2.  SCOPE 

This  section  (Volume  II)  of  the  report  describes  the  workings  of  the 
second  program,  PHASE2.  A program  overview  Is  given,  followed  by  a thorough 
description  of  the  logic  Involved  In  map  processing.  A skeleton  of  the  logic 
flow  Is  provided.  Input  and  output  files  are  described.  Program  require- 
ments and  restrictions,  error  messages  and  error  handling  techniques,  defi- 
nitions of  Import  variables,  and  an  estimate  of  running  time  are  also  pre- 
sented . 
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SECTION  II 
PROGRAM  OVERVIEW 


Determining  the  minimum  number  of  collection  vehicles  needed  to  service 
a base  is  fairly  simple;  the  task  of  assigning  collection  schedules  in  such 
a way  that  each  vehicle  is  used  to  capacity,  but  not  overfilled,  while  travel 
time  and  distance  are  kept  close  to  the  minimum,  is  more  difficult. 

Program  PHASE2  serves  two  purposes:  it  assigns  street  segments  to  sec- 
tions (a  section  is  a set  of  streets  to  be  serviced  by  one  collection  vehicle), 
and  it  plots  the  results  of  the  sectioning.  The  sectioning  groups  the  street 
segments  into  closely  connected,  reasonably  convex  sets.  The  size  of  each 
section  Is  determined  by  the  capacity  of  the  refuse-collection  vehicle  that 
will  service  It.  Therefore,  choosing  the  streets  in  a section  so  that  each 
section  is  compact  is  the  main  effort  in  PHASE2. 

Two  types  of  data  are  used  as  input  to  PHASE2.  Data  describing  the  nodes 
and  segments  are  read  from  files  TAPEll  and  TAPE9.  Card  input  is  used  to  spe- 
cify the  problem  title,  the  vehicle  capacities  and  numbers,  the  time  limits, 
the  base  segment  for  the  first  section,  and  the  map  bounds. 

The  program  consists  of  a main  program,  PHASE2,  and  11  subroutines. 

PHASE2  reads  the  data  cards,  the  segment  data  from  file  TAPE9,  and  the  node 
data  from  file  TAPEll.  Refuse-quantity  Information  included  with  the  segment 
data  and  vehicle-capacity  data  from  the  data  cards  are  used  to  determine  the 
number  of  vehicles  required  to  collect  all  of  the  refuse. 

PHASE2  calls  subroutine  BUILD  to  build  a near-neighbor  table.  The  table 
indicates,  for  each  street  segment,  the  60  other  segments  closest  to  it.  To 
build  the  table,  subroutine  BUILD  computes  the  distances  from  each  segment  to 
each  other  segment.  The  60  shortest  distances  and  the  corresponding  segment 
numbers  are  found  using  an  in-core  tree-sort  algorithm  in  subroutines  SORTK 
and  SIFTUP. 


PHASE2  calls  subroutine  SECTION  to  assign  the  segments  to  sections  (cot 
responding  to  collection  trips).  Segments  are  selected  for  addition  to  a 
section  on  the  basis  of  the  number  of  near  neighbors  they  share  with  another 
segment,  called  a base  segment,  already  in  the  section.  The  first  base  seg- 
ment is  specified  by  the  user.  Subsequent  base  segments  are  selected  as  the 
sections  are  built.  Segments  are  added  to  a section  as  long  as  the  vehicle 
capacity  is  not  exceeded. 

After  each  section  is  filled,  subroutine  SECTION  checks  to  see  whether 
all  remaining  unassigned  segments  will  fit  into  a single,  last  section  If 
not,  selection  continues  on  the  basis  of  the  shared  near-neighbor  criterion. 

As  sections  are  completed,  the  segment  numbers  are  written  to  file  TAPE4, 
and  statistics  on  vehicle  time  and  capacity  are  accumulated. 

After  the  sectioning  has  been  completed,  PHASE2  calls  subroutine  PLOTS 
to  initialize  the  plotting  package.  Subroutine  MAPPLT  is  called  to  plot  maps 
indicating  the  section  assignment  of  each  segment.  MAPPLT  uses  subroutine 
SHAPCOM  to  set  shape  parameters  for  the  segments  and  subroutine  COORD  to  gen- 
erate coordinates  of  points  on  each  segment.  Subroutine  NUMBER  is  called  to 
append  section  numbers  to  the  segments.  Program  PHASE2  terminates  after  a 
trip  data  summary  is  written  to  file  TAPE1. 

The  flow  of  control  from  one  subprogram  to  another  is  shown  in  Figure  1. 
Within  each  subprogram,  only  the  first  call  to  each  other  subprogram  is  shown. 
(Three  of  the  subroutines  shown  in  Figure  1--PL0TS,  PLOT,  and  SYMBOL--  are 
subroutines  from  the  basic  Calcomp  software  package  and  are  not  included  in 
ihe  description  of  program  PHASE2.) 
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SECTION  III 
PROGRAM  LOGIC 


The  logic  for  program  PHASE2  is  described  from  three  viewpoints.  The 
first  description  is  task  oriented.  The  second  is  data-storage  oriented  and 
includes  discussions  of  the  preparation  of  data  for  use  by  subsequent  programs, 
the  use  of  input  data,  and  the  data  structures  used  in  PHASE2.  The  third  view 
describes  each  subroutine  in  terms  of  its  purpose  and  the  manipulations  per- 
formed within  it. 

1 . PROGRAM  TASKS 

Three  tasks  are  accomplished  by  program  PHASE2:  (1)  the  number  of  trips 
is  determined  on  the  basis  of  one  of  two  options  available  to  the  user;  (2) 
street  segments  are  assigned  (in  accordance  with  vehicle-capacity  restrictions) 
to  sections,  each  section  corresponding  to  a collection-vehicle  trip,  on  the 
basis  of  the  number  of  near  neighbors  each  segment  shares  with  some  segment 
already  in  the  section;  (3)  finally,  one  or  more  maps  are  plotted  showing  the 
section  assignments  of  the  segments. 

Program  execution  begins  in  the  main  program,  PHASE2.  The  problem  title 
is  read  from  the  first  data  card.  The  number  of  vehicles  and  their  capacities, 
the  time  limits,  and  the  number  of  the  segment  that  is  to  be  the  first  base 
segment  are  read  from  the  next  two  data  cards.  The  number  of  vehicles  speci- 
fied on  the  data  cards  determines  the  method  used  to  generate  the  number  of 
trips.  If  enough  vehicles  are  specified  to  collect  all  the  refuse  in  the  col- 
lection region,  then  that  number  of  vehicles  is  used,  even  if  it  is  not  the 
minimum.  If  fewer  vehicles  are  specified  than  are  needed  to  service  the  en- 
tire region,  the  program  will  assign  vehicles  in  the  order  given  on  the  data 
cards  until  the  minimum  number  needed  to  collect  all  the  refuse  is  obtained. 

Segment  data  are  read  from  file  TAPE9,  and  refuse-quantity  and  node 
data  are  read  from  file  TAPE  11.  The  input  data  and  the  number  of  vehicles 
that  will  be  needed  are  printed.  Subroutine  BUILD  is  called  to  construct 
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.*  net>  -neighbor  table,  fo1  ea  i segment,  BUI1D  computes  the  distances  4 
each  other  segment.  The  nun be> s of  the  other  segments  are  masked  into  th 
low-order  12  bits  of  the  distance.  The  distances  are  tree-sorted  by  sub- 

■ out  ices  SORTK  and  S1FTUP.  Tl'e  tee  )'t  orders  only  t tie  spec* Lied  r be-  of 
i te  ' The  60  shortest  distances  ire  obtained  from  subroutine  SORTK,  and 

the  segment  numbers  are  retrieved  from  these  distances.  Individual  bits 
oi responding  to  the  segment  numbers  are  set  to  1 in  an  array  called  a near- 
neighbor  list.  The  near-neighbor  list  and  information  describing  the  origi- 

■ r segment  are  written  to  disk  This  procedure  is  repeated  for  each  segment 
in  the  map  PHASE2  then  calls  subroutine  SECTION  to  assign  segments  to 
sections. 

Subroutine  SECTION  chooses  segments  to  be  added  to  a section  by  deter 

■ icing  which  segments  share  the  most  near  neighbors  with  a base  segment  in 
the  section.  Segments  are  added  to  the  section  as  long  as  the  vehicle's 
capacity  and  time  limit  are  not  exceeded.  A section  is  complete  when  a cer- 

airi  minimum  load  has  been  achieved  or  when  each  remaining  segment  woe  causa 
the  vehicle's  capacity  or  time  limit  to  be  exceeded.  The  minimum- load  cri- 
*0*v  rie?  te  make  the  ctmtJ-aiive  load  at  the  completion  of  a section  eeeal 
to  that  fraction  of  the  total  refuse  corresponding  to  the  ratio  of  vehicle 
capacity  used  to  total  vehicle  capacity  available.  If  one  or  more  sections  j 
. c closed  out  because  each  of  the  remaining  segments  would  cause  the  vehicle 
capacity  to  be  exceeded,  a situation  may  occur  wherein  the  determined  nruimum 
number  of  vehicles  will  be  inadequate  to  collect  all  of  the  refuse.  Tn  tiiis 
c.’so  additional  vehicles  are  assigned  in  the  order  in  which  the  vehicles  have 
Lee  specified,  and  a message  is  printed.  As  each  section  is  completed,  the 
■ et  lining  refuse  quantity  is  determined;  if  the  remaining  segments  can  be 
assigned  to  one  vehicle,  the  shared  near-neighbor  testing  is  discontinued 
a.-:  . : of  the  remaining  segments  are  assigned  to  the  last  section. 

i 

I 

After  subroutine  SECTION  has  completed  the  sectioning,  PHASE2  prints  a 
summary  of  the  loads  and  times  required  by  the  vehicles.  A list  of  the 
numbers  of  the  segments  in  each  section  is  also  printed. 

i 

i 

Subroutine  PLOTS  is  called  to  initialize  the  plotting  package.  Subrou- 
tine PLOT  is  called  to  place  a 3-inch  border  at  the  bottom  of  the  plot. 
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PHASE2  then  reads  map  bounds  from  the  remaining  data  cards.  If  no  map- 
bounds  cards  are  found,  defaults  are  set  up  for  a 30-  by  30-inch  map. 

Subroutine  MAPPLT  is  called  once  per  output  map.  MAPPLT  examines  sequen- 
tially the  original  segment  data,  skipping  segments  that  are  outside  the  map 
bounds  and  drawing  segments  that  lie  at  least  partially  within  bounds.  Before 
each  segment  is  drawn,  subroutine  SHAPCOM  is  called  to  set  up  parameters  used 
to  determine  the  position  on  the  segment  of  points  at  a given  distance  from 
the  start  of  that  segment.  The  actual  coordinates  of  the  points  on  the  segment 
are  returned  by  subroutine  COORD. 

Four  different  computations  are  used  by  subroutines  SHAPCOM  and  COORD  to 
produce  the  coordinates  of  points  on  a segment.  The  simplest  computation  Is 
performed  for  straight  segments  and  involves  a linear  interpolation  between  the 
initial  and  final  nodes.  Another  calculation  processes  both  circular-arc  and 
S-curve  segments;  an  S-curve  is  treated  as  two  consecutive  circular  arcs.  In 
calculating  the  coordinates  of  a point  on  a rectangular  segment,  the  Jope  com- 
ponents of  the  first  side  are  determined;  appropriate  multiples  of  these  com- 
ponents are  then  added  to  the  starting  or  ending  node's  coordinates.  The  coor- 
dinates of  a point  on  an  angle  segment  are  found  by  linear  interpolation  be- 
tween one  end  of  the  angle  and  the  vertex.  (A  full  description  of  the  geometry, 
as  well  as  relevant  calculations,  are  given  in  Section  III  of  Volume  I of 
this  report.  The  scale  ratio  SCR  is  replaced  by  1.0  in  program  PHASE2.) 

2.  DATA  STORAGE 

Program  PHASE2  obtains  data  from  three  sources:  card  input,  file  TAPE9, 
and  file  TAPE11.  Two  files,  TAPE1  and  TAPE4,  are  generated  by  program  PHASE2 
and  are  saved  on  disk  for  use  by  program  PHASE3.  Files  TAPE2,  TAPE3,  TAPE7, 
and  TAPE10  are  used  as  scratch  files. 

The  card  data,  the  data  from  TAPE9,  and  the  data  from  TAPE11  are  read  at 
the  beginning  of  PHASE2.  The  segment  data  read  from  TAPE9  are  stored  In  arrays 
in  blank  COMMON.  The  street  number  and  the  number  of  ways  of  travel  on  TAPE9 
are  not  retained  in  core.  The  node  data  from  TAPE 11  are  stored  in  arrays  in 


lain* led  COMMON  block  NOOATA.  It  tin*  number  of  vehicles  from  the  second  data 
catil  is  zero  for  any  vehicle,  it  is  reset  to  1 in  the  loop  on  statement  ?‘j. 

'be  numbers  of  vehicles,  their  capacities,  and  their  time  limits  are  moved  to 
the  I ROCKS  array  in  the  loop  through  statement.  SO.  liu*  amount  of  total  refuse 
*s  also  accumulated  in  this  loop.  When  subroutine  HUILl)  is  called  to  build  the 
"eat  neighbor  table,  segment  numbers  in  array  1STPR  and  street  segment  midpoints 
in  arrays  X,  Y,  XT,  and  YT  are  sent  through  the  argument  list.  Subroutine  lUMl.D 
writes  to  l AH'/  the  segment  number,  refuse  quantity,  travel  and  collection  time, 
number  of  houses,  and  .’6  words  of  near  neighbor  information  for  each  segment. 

When  PHASE  2 calls  subroutine  SECTION,  the  vehicle  data  in  array  1 RUCKS, 
ditch  are  grouped  by  vehicle  capacity,  are  expanded  into  array  I RUCK  so  that 
each  line  of  array  TRUCK  corresponds  to  a single  vehicle.  As  the  first  base 
• eqment  is  sought  in  the  loop  through  statement  1/,  segments  other  than  t tie 
base  segment  are  written  to  lAI'l  1.  1 AIM  1 will  contain  unassigned  segments 

and  then  near  neighbor  lists.  As  each  segment  is  considered  tor  addition 
to  a section,  its  segment  and  neighbor  data  are  written  to  lAI't?  if  it  is  not 
nided  to  the  section.  As  a segment  is  assigned  to  a section,  its  segment  and 
neighbor  data  are  written  to  tile  1APE.T.  When  a section  is  completed,  the 
segment  numbers  are  read  from  file  IAPt.1  and  written  to  file  TAPE4,  The  un- 
assigned segments  on  file  TAPI  1 are  recopied  to  file  TAPE?.  When  all  of  the 
seqments  have  been  assigned  to  sections,  control  returns  to  program  PHAM?. 

\t  ttii'  end  of  program  PHASE?,  the  number  of  segments  and  the  number  of  sec 
tions  are  written  to  IAI'11.  as  are  pointers  to  the  first  and  last  segment 
numbei  s on  tilt'  IAPM  and  the'  vehicle  rapacity  tor  each  section. 


PURPOSl  AND  PFKt  ORMANCf 

In  this  section  the  simpler  subroutines  are  described  first  so  their 
workings  will  be  clear  when  they  are  mentioned  again  in  the  descriptions  of 
the  more  complicated  subroutines  and,  finally,  of  the  main  program.  Logic 
flowcharts  are  given  m Appendix  A.  Complete  program  listings  are  provided 
in  Appendix  It.  In  Appendix  C,  the  more  important  variables  mentioned  in  the 
following  descriptions  are  defined  in  terms  of  their  specific  meaning  for 
each  subroutine. 
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a.  Function  KOUNT 


The  purpose  of  function  KOUNT  is  to  count  the  1 bits  in  a 60-bit  word. 
The  function  is  written  in  the  COMPASS  assembler  language  for  the  CDC  6600. 

Function  KOUNT  has  one  argument.  The  argument  is  a bit  pattern 
with  bits  set  to  1 where  two  segments  share  a near  neighbor.  The  value  re- 
turned is  the  number  of  1 bits  in  the  argument.  The  sequence  of  instructions 
generated  by  the  compiler  where  KOUNT  is  called  includes  setting  register  A1 
to  the  address  of  the  argument  list.  In  KOUNT,  the  first  SA1  Instruction 
causes  the  XI  register  to  receive  the  address  of  the  first  argument.  The 
second  SA1  instruction  causes  the  value  of  the  argument  to  be  placed  in  reg- 
ister XI.  The  CX6  instruction  counts  the  1 bits  in  register  XI  and  places 
the  result  in  register  X6.  Control  then  returns  to  the  calling  program. 

b.  Subroutine  SHLSRT 

Subroutine  SHLSRT  sorts  one  array  into  decreasing  order  and  carries 
a second  array  along  during  the  sorting.  The  algorithm  used  is  Shell's  sort- 
ing algorithm. 

Subroutine  SHLSRT  has  three  arguments.  The  first  argument  is  the 
array  to  be  sorted.  The  second  argument  is  an  array  that  is  paired  with  the 
array  to  be  sorted  and  is  rearranged  as  the  first  array  is  sorted.  The  third 
argument  is  the  number  of  words  to  be  sorted. 

The  statements  up  to  statement  60  arrange  array  X in  Increasing 
order.  The  numbers  are  sorted  by  a procedure  in  which  pairs  of  numbers  are 
compared  and  Interchanged  If  necessary  to  put  the  smaller  number  closer  to 
the  beginning  of  the  array.  The  separation  of  the  numbers  compared  Is  approx- 
imately one-half  the  number  of  entries  in  the  array;  this  spacing  is  halved  in 
subsequent  passes  through  the  array.  When  two  numbers  are  Interchanged,  the 
pointers  are  moved  up  so  that  the  smaller  number  is  compared  to  a number  far- 
ther up  in  the  array.  The  spacing  (N)  is  set  initially  to  one-half  the  number 
of  words.  The  number  of  comparisons  (K)  to  be  performed  in  the  loop  through 
statement  50  is  computed  as  the  total  number  of  words  less  N. 
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lire  loop  through  statement  SO  uses  index  1 as  one  of  the  pointers, 
lhis  pointer  is  sorted  in  variable  J.  The  other  pointer,  L,  is  set  equal  to 
l*N.  I he  values  of  the  array  to  he  sorted  and  the  array  to  be  carried  alomi 
are  saved  as  XT  and  A1 . The  values  of  X at  the  locations  indicated  by  the 
pointers  are  compared;  If  they  are  in  order,  control  transfers  to  statement  40. 
If  not,  the  larger  value  Is  stored  closer  to  the  end  of  the  array.  The  pointers 
are  both  moved  up  by  N;  if  the  smaller  valued  pointer  is  a valid  subscript,  con- 
trol transfers  to  statement  20,  where  another  comparison  is  performed. 

At  statement  40  the  saved  values  are  stored  in  the  appropriate  place 
in  the  arrays.  When  the  loop  through  statement  50  is  completed,  if  the  spac- 
ing is  equal  to  1,  the  sort  is  complete  and  control  transfers  to  statement  60. 
Otherwise,  t tie  spacing  is  halved  and  control  transfers  to  statement  10. 

I he  loop  through  statement  70  rearranges  the  arrays  so  that  the  X- 
arrav  is  in  decreasing  order.  Control  returns  to  the  calling  program. 

c.  Subroutine  SIFTUP 

Subroutine  SIFTUP  orders  each  subtree  in  a binary  tree  from  a given 
subroot  up  to  the  root  so  that  each  subroot  is  smaller  than  either  of  its 
branches.  Subroutine  SIFTUP  has  four  arguments.  The  first  is  the  subscript 
of  the  subroot  at  which  sorting  starts.  The  second  is  the  number  of  items 
in  the  array  to  be  sorted.  The  third  is  the  array  to  be  sorted.  The  fourth 
is  the  dimension  of  tire  array  to  be  sorted. 

Variable  I is  set  equal  to  the  subscript  of  the  subroot  where  sort 
' rig  will  start.  Variable  I will  continue  to  point  to  a subroot  throughout 
the  sorting.  I he  value  at  TRFF(l)  is  saved  in  variable  COPY. 


At  statement  10,  pointer  J is  set  equal  to  the  subscript  of  the  left 
branch.  It  .1  points  to  or  beyond  the  last  item  in  the  tree,  control  transfers 
to  statement  6.  Otherwise,  the  left  and  right  branches  are  compared.  If  the 
right  branch  Is  smaller,  J will  be  incremented  so  that  it  points  to  the  right 
br  am  h . 


At  statement  4 the  smaller  branch  Is  compared  to  the  root.  If  the 


root  is  smaller,  control  transfers  to  statement  6.  Otherwise,  control  resumes 
at  statement  5,  and  the  branch  value  is  stored  in  the  root  position.  The 
branch  from  which  the  smaller  number  came  now  becomes  the  root  in  another  Iter- 
ation. Control  transfers  back  to  statement  10. 

At  statement  6,  the  value  of  the  root  saved  in  variable  COPY  is 
stored  in  the  appropriate  place  in  array  TREE.  Control  returns  to  the  call- 
ing program. 

d.  Subroutine  SORTK 

Subroutine  SORTK  returns  the  KN  smallest  numbers  In  array  TREE.  A 
tree-sort  algorithm  is  used.  The  array  to  be  sorted  is  treated  as  a binary 
tree  and  is  partially  ordered,  in  such  a manner  that  each  subroot  Is  smaller 
than  its  branches,  by  calls  to  subroutine  SIFTUP. 

V 

Subroutine  SORTK  has  four  arguments.  The  first  (N)  is  the  number  of 
items  in  the  array.  The  second  (KN)  gives  the  number  of  items  to  be  returned. 

The  third  (TREE)  Is  the  array  to  be  sorted.  The  fourth  is  the  dimension  of 
the  array. 

Subroutine  SORTK  begins  by  comparing  the  first  and  last  numbers  in 
the  array  to  be  sorted.  If  the  last  item  is  smaller  than  the  first,  the  two 
are  interchanged  to  prevent  the  smaller  number  from  being  trapped  as  the  last 
entry  in  the  array  when  the  number  of  entries  is  even  and  the  last  entry  is 
the  smallest.  Variable  K is  set  to  one-half  the  number  of  items  in  the  tree. 

The  loop  through  statement  10  calls  SIFTUP;  the  first  argument  starts 
at  the  middle  of  the  tree  and  works  back  to  the  second  element  In  the  array. 

When  this  loop  Is  complete,  the  branches  in  the  tree  are  smaller  than  any  sub- 
roots except  the  root  of  the  entire  tree. 

The  loop  through  statement  11  causes  SIFTUP  to  move  the  smallest  num- 
ber to  the  root  of  the  tree.  This  number  Is  then  exchanged  for  the  last  Item 
in  the  tree.  The  loop  through  statement  11  is  used  once  for  each  number  to  be 
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returned.  When  this  loop  is  complete,  the  KN  smallest  numbers  will  be  at  the 
end  of  array  TREE,  and  the  smallest  number  will  be  last.  Control  then  returns 
to  the  calling  program. 

e.  Function  IFIND 

Function  IFIND  uses  a binary  search  to  locate  a given  number  in  an 
art  ay;  the  subscript  corresponding  to  the  location  of  the  number  is  assigned 
as  the  value  of  IFIND.  If  the  number  is  not  found,  the  function  sets  the 
value  of  IFIND  equal  to  the  negative  of  the  subscript  at  which  the  number,  to 
be  in  numerical  order,  should  be  inserted.  (The  array  is  assumed  to  be  in  in- 
creasing order.) 

The  comment  cards  at  the  beginning  of  function  IFIND  list  the  latest 
changes  to  the  function  and  state  the  function's  purpose. 

Argument  NUM  is  the  number  that  is  sought  in  array  IARRAY . The 
length  of  array  IARRAY  is  given  by  argument  LEN.  Function  IFIND  begins  by  check- 
ing that  LEN  ^ 0.  If  LEN  < 0,  the  function  assigns  a value  of  -1  to  IFIND.  This 
value  indicates  that  the  number  sought  is  not  in  the  array  and  would  be  stored 
as  the  first  entry  in  the  array.  The  binary  search  uses  variables  II,  IP,  and 
IF  as  pointers.  II  is  the  subscript  of  the  front  of  the  region  being  searched, 
IP  is  the  subscript  of  the  item  being  compared  to  the  number  sought,  and  IF  is 
the  subscript  of  the  last  item  in  the  region  being  searched.  Variable  II  is 
initially  set  to  1 at  statement  5,  and  variable  IF  is  set  to  the  end  of  the 
array  in  the  next  statement.  The  pointer,  IP,  is  the  subscript  about  midway 
between  II  and  IF. 

The  computation  of  IP  occurs  at  statement  10.  The  statement  follow- 
ing statement  10  compares  the  number  being  sought,  NUM,  to  the  data  at 
IARRAY (IP).  If  NUM  < IARRAY(IP),  control  transfers  to  statement  20,  indica- 
ting that  the  number  is  in  the  front  one-half  of  the  region  being  searched;  at 
statement  20  the  final  pointer  is  moved  to  the  subscript  preceding  the  point 
just  searched.  If  NUM  > IARRAY (IP),  control  transfers  to  statement  30,  indi- 
cating that  the  number  being  sought  follows  the  subscript  just  inspected.  At 
statement  30  the  initial  pointer,  II,  is  set  to  the  present  pointer,  IP,  plus  1. 
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If  the  number  sought  Is  found  at  IARRAY(IP),  control  transfers  to  statement  50, 
where  IFIND  Is  set  equal  to  the  current  pointer  and  control  returns  to  the  call- 
ing program.  Where  NUM  Is  unequal  to  IARRAY(IP),  control  resumes  at  statement 
40  after  the  initial  or  final  pointers  are  moved.  At  statement  40  the  final 
pointer  is  compared  to  the  initial  pointer;  if  IF  > II,  control  is  transferred 
to  statement  10. 

At  statement  10  the  search  is  resumed  on  the  appropriate  one-half  of 
the  region  examined  previously.  If  the  final  pointer  becomes  less  than  the 
initial  pointer,  the  number  sought  is  not  in  the  array.  In  this  case,  control 
resumes  following  statement  40,  and  the  value  of  IFIND  is  set  to  the  negative 
of  the  current  pointer.  If  the  number  at  the  current  pointer  is  less  than  the 
number  being  sought,  IFIND  is  set  to  -(IP  + 1 ) so  the  number  can  be  inserted 
in  the  appropriate  place.  Control  then  returns  to  the  calling  program. 

f.  Subroutine  NUMBER 

Subroutine  NUMBER  appends  numbers  to  plotted  output.  Its  purpose  is 
almost  identical  to  that  of  the  standard  Calcomp  number  routine,  the  primary 
difference  being  that  the  last  argument  in  subroutine  NUMBER  gives  an  alpha- 
numeric format  rather  than  an  integer  format  code. 

Subroutine  NUMBER  has  six  arguments.  The  first  two  give  the  coor- 
dinates, in  plotter  inches,  of  the  lower  left  corner  of  the  field.  The  third 
gives  the  height,  in  inches,  of  the  digits.  The  fourth  is  the  number  to  be 
plotted.  The  fifth  is  the  angle  at  which  the  number  is  to  be  plotted,  measured 
in  degrees  counterclockwise  from  the  horizontal.  The  last  argument  is  an  alpha- 
numeric format  up  to  10  characters  long,  which  describes  the  appearance  of  the 
plotted  number. 

Array  TEXT  is  used  to  hold  the  character  representation  of  the  number. 
Up  to  30  characters  are  allowed.  The  first  executable  FORTRAN  statement  sets 
this  array  to  three  words  of  blanks.  The  second  statement  moves  the  format  in- 
to the  second  word  of  array  FORM.  The  first  and  third  words  of  this  array  have 
been  preset  to  a left  and  a right  parenthesis  by  a DATA  statement.  The  ENCODE 
statement  converts  the  number  from  binary  form  in  variable  NUM  to  character 
form  in  array  TEXT,  according  to  format  FORM. 
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A character  count,  variable  NC,  is  set  to  30.  The  loop  through  state- 
ment 10  searches  for  the  last  non-blank  character  in  array  TEXT.  Each  time  a 
blank  is  found,  starting  at  the  end  of  the  TEXT  array,  the  character  count  (NC) 
is  decremented  by  1.  When  a non-blank  character  is  encountered,  control  trans- 
fers to  statement  20.  Statement  20  calls  the  standard  SYMBOL  subroutine  to 
plot  the  character  representation  of  the  number.  Control  then  returns  to  the 
calling  program. 

g.  Subioutine  SHAPCOM 

Subroutine  SHAPCOM  sets  up  parameters  in  COMMON  block  COPARM  that  de- 
scribe the  geometrical  properties  of  a segment.  These  parameters  are  used  by 
subroutine  C00R0  to  produce  the  coordinates  of  points  on  a segment. 

Subroutine  SHAPCOM  has  two  arguments.  Argument  TOTLEN  gives  the 
total  length  of  the  segment,  in  miles.  Argument  AVMD  gives  the  number  of  miles 
per  map  coordinate  unit  (MCU)  on  the  first  map  input  to  program  RCINPT.  The 

values  of  the  arguments  are  sent  to  subroutine  SHAPCOM,  and  all  output  values 

from  SHAPCOM  are  placed  in  COMMON  block  COPARM. 

In  COMMON  block  COPARM,  variable  SF  indicates  the  shape  of  the  seg- 
ment. XN I and  XNF  are  the  x-coordinates  of  the  initial  and  final  nodes  of  the 
segment.  YN1  and  YNF  are  the  y-coordinates  of  these  nodes.  SX  and  SY  are  the 
slope,  in  MCU  per  mile,  in  the  x and  y directions.  RPR  is  the  »'eciprocal  of 

the  radius  of  curvature  for  circular  segments  and  the  circular  portions  of  S- 

curves.  Cll  and  Cl 2 are  the  position  differences,  in  MCU,  of  the  starting 
point  and  center  of  a circular  arc  or  of  the  first  one-half  of  an  S-curve. 

XCTR  and  YCTR  are  the  center  coordinates,  in  MCU,  for  a circular  arc  or  one- 
half  of  an  S-curve.  BR1  is  the  distance  in  miles  from  the  start  of  a segment 
to  some  particular  point  on  that  segment.  It  is  not  used  for  straight  segments. 
For  circular  segments,  BR1  is  the  total  perimeter.  For  an  S-curve,  BR1  is 
the  perimeter  to  the  midpoint  of  the  S-curve.  For  a rectangular  segment,  BR1 
is  the  distance  to  the  first  bend  in  the  rectangle.  For  an  angle,  BR1  is  the 
distance  to  the  vertex.  BR2  is  defined  only  for  rectangular  segments  and 
angles.  For  a rectangular  segment,  it  is  the  perimeter  in  miles  from  the  start 
of  the  segment  to  the  second  bend.  For  an  angle,  BR2  is  the  length  of  the  sec- 
ond side.  SGN  is  -1  for  shapes  involving  the  L (left)  prefix,  and  +1  otherwise. 
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Subroutine  SHAPCOM  begins  execution  by  assuming  that  the  shape  code 
indicates  a straight  line.  Break  indicators  BR1  and  BR2  are  set  to  0.  DX  and 
DY , the  x-  and  y-components  of  the  vector  from  the  initial  to  the  final  node 
on  the  segment,  are  computed.  The  x-  and  y-components  of  the  slope  of  the  vec- 
tor, measured  in  MCU  per  mile,  are  computed  and  stored  in  SX  and  SY.  The  shape 
code  is  tested;  if  the  segment  proves  to  be  a straight  line  or  is  not  to  be 
plotted,  the  subroutine  returns  control  to  the  calling  program.  For  any  other 
shape  code,  execution  continues.  The  angle  of  the  vector  from  the  starting  to 
the  stopping  node  is  computed  as  variable  THETA.  The  distance  from  the  start- 
ing to  the  stopping  point,  D,  is  computed  in  miles.  If  the  shape  code  indicates 
a shape  other  than  circular  or  S-curve,  control  transfers  to  statement  60. 


At  statement  45  the  coordinates  of  the  final  node  are  stored  in  vari- 
ables XE  and  YE.  The  first  break,  BR1 , is  set  to  the  total  length  of  the  seg- 
ment. Variable  DD  is  set  to  the  straight-line  distance  from  the  starting  to 
the  stopping  node.  If  the  shape  code  indicates  a circular  segment,  control 
transfers  to  statement  50.  If  not,  variables  XE  and  YE  are  reset  to  the  coor- 
dinates of  the  midpoint  of  the  S-curve.  Break  indicator  BR1  is  reset  to  the 
perimeter  length  from  the  starting  point  to  the  center  of  the  S-curve.  Variable 
DD  is  set  to  one-half  the  distance  from  the  starting  to  the  stopping  point. 

At  statement  50,  SGN  is  set  to  1 . If  the  shape  code  indicates  a 
left  circle  or  left  S-curve,  SGN  is  reset  to  -1.  Variable  V is  set  equal  to 
1-D/T0TLEN.  VS  is  the  square  of  V.  The  reciprocal  of  the  radius  of  curva- 
ture of  the  circle  or  the  circular  portion  of  the  S-curve  is  evaluated  using 
a polynominal  approximation  to  the  solution  from  a transcendental  equation 
containing  the  reciprocal  of  the  radius  of  curvature.  The  approximate  radius 
of  curvature,  RPR,  is  improved  by  a series  of  linear  interpolations  if  the 
value  for  RPR  causes  an  error  greater  than  0.00001  in  the  transcendental 
equation 

BR1*RPR  DD*RPR 

sin  ^ — — 2^ — 

When  RPR  is  within  the  desired  accuracy,  control  resumes  at  state- 
ment 51.  The  radius  of  curvature,  R,  is  computed.  A temporary  variable, 

ARG,  is  evaluated.  The  height  of  the  center  of  the  circle  from  the  line 
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connecting  the  starting  and  stopping  points,  H,  is  set  to  0.  If  variable  ARG 
is  greater  than  0,  H is  recomputed.  The  distance  to  the  first  break,  BR1 , is 
tested  to  determine  whether  the  circular  arc  is  greater  than  one-halt  a circle 
If  so,  the  sign  of  the  height  is  changed.  The  x-  and  y-coordinates  of  the 
center  of  the  circle  are  computed.  The  components  of  the  vector  from  the  cen- 
ter to  the  starting  point,  Cll  and  C 1 2 , are  computed.  All  variables  needed  to 
compute  points  on  the  S-curve  or  circle  are  now  available,  so  control  returns 
to  the  callimi  program. 

Processing  continues  at  statement  60  for  the  remaining  shape  codes. 
At  statement  60,  the  shape  code  is  tested;  if  neither  a right  nor  a left  rec- 
tangle is  indicated,  control  transfers  to  statement  80.  Otherwise,  for  a rec- 
tangular segment,  the  distance  from  the  start  to  the  first  bend,  BK 1 , is  com- 
puted. If  this  distance  is  greater  than  0.05  of  the  total  length,  control 
transfers  to  statement  70.  Otherwise,  the  rectangle  is  assumed  to  be  so  shal- 
low that  a straight-line  approximation  is  adequate,  and  the  shape  code  is  set 
to  0.  Control  then  returns  to  the  calling  program. 

At  statement  70  the  perimeter  to  the  second  bend  in  the  rectangle, 
is  computed.  SX  and  SY , the  x-  and  y- components  of  the  slope  of  the 
vector  from  starting  point  to  stopping  point,  are  computed,  and  control  re- 
turns to  the  calling  program. 

The  only  segments  that  reach  statement  80  are  the  angles.  The  sign 
of  the  shape  code  is  retrieved  in  variable  SON,  and  the  distance  from  the 
starting  node  to  the  vertex  of  the  angle  is  retrieved  as  the  magnitude  of  the 
shape  code  and  is  stored  in  variable  BR 1 . The  length  of  the  second  leg  of  the 
angle  is  computed  and  saved  in  variable  BR2.  If  the  angle  is  incorrectly  spe- 
lfied  so  that  it  is  actually  a straight  segment,  a round-off  error  may  occur 
n the  computation  of  ARG,  the  square  of  the  distance  from  the  vertex  to  the 
line  connecting  the  endpoints.  If  ARG  is  zero  or  negative,  control  transfers 
to  statement  100.  Otherwise,  the  x-  and  y-coordinates  of  the  vertex  are  com- 
puted. and  control  returns  to  the  calling  program. 

At  statement  100,  the  shape  code  and  break  indicators  are  set  to 
O,  indicating  a straight  segment.  Control  returns  to  the  calling  program. 
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h.  Subroutine  COORD 


■ 


Subroutine  COORD  Is  given  a distance,  In  miles,  from  the  beginning  of 
a segment  and  returns  the  coordinates  In  MCU.  Parameters  describing  the  segment 
to  be  processed  have  been  sorted  In  COMMON  block  COPARM  by  subroutine  SHAPCOM 
before  COORD  Is  called.  Argument  CUMLEN  Is  the  cumulative  length  along  the 
string.  In  miles;  arguments  XX  and  YY  are  the  coordinates  returned  for  a point 
CUMLEN  miles  from  the  start  of  the  segment. 

The  first  statement  of  COORD  sets  S equal  to  the  cumulative  length. 

If  the  shape  code  Is  nonzero,  control  transfers  to  statement  10.  The  coordin- 
ates of  the  point  on  a straight-line  segment  are  computed  and  returned  In  vari- 
ables XX  and  YY.  Control  returns  to  the  calling  program. 

At  statement  10  control  transfers  to  statement  30  If  the  shape  code 
Indicates  other  than  a circular  or  S-curve  segment.  For  circular  and  S-curve 
segments,  the  reciprocal  of  the  radius  of  curvature  Is  stored  In  RIP.  The  co- 
ordinates of  the  center  of  the  circular  portion  are  stored  In  XC  and  YC.  The 
components  of  the  vector  from  the  center  of  the  circle  to  the  Initial  node  are 
stored  In  Cl  and  C2.  If  the  point  on  the  segment  Is  less  than  or  equal  to 
0.999  of  the  first  break  distance,  or  If  the  shape  code  Indicates  a circular 
segment,  control  transfers  to  statement  20.  The  statements  following  this 
test  change  parameters  to  generate  coordinates  for  the  second  circular  portion 
of  an  S-curve.  The  sign  of  the  reciprocal  of  the  radius  of  curvature  Is  re- 
versed. The  cumulative  distance,  S,  Is  set  to  the  distance  from  the  midpoint 
of  the  S-curve.  The  coordinates  of  the  center  of  the  second  circular  portion, 

XC  and  YC,  are  computed.  Variables  Cl  and  C2  are  recomputed  for  the  new  center. 

At  statement  20  the  sine  and  cosine  of  the  angle  subtended  by  the 
perimeter  corresponding  to  S are  computed.  The  coordinates  of  the  point,  XX 
and  YY,  are  computed,  and  control  returns  to  the  calling  program. 

At  statement  30,  control  transfers  to  statement  60  If  the  shape  code 
Indicates  that  the  segment  Is  not  a rectangle.  Otherwise,  variable  SON  Is  set 
to  1.  If  the  shape  code  Indicates  a left  rectangle,  SON  Is  reset  to  -1.  If  S, 
the  distance  along  the  rectangle,  Is  greater  than  1.05  times  the  first  side's 
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length,  control  transfers  to  statement  40.  If  S is  greater  than  0.95  times 
the  length  of  the  first  leg,  S is  set  to  the  length  of  the  first  leg.  The  x- 

and  y-coordlnates  of  the  point  on  the  first  leg  are  computed  by  linear  inter- 

polation, and  control  returns  to  the  calling  program. 

At  statement  40,  S Is  tested  to  see  whether  It  falls  on  the  second 
leg  of  the  rectangle.  If  S Is  greater  than  1.05  times  BR2,  the  length  of  the 
second  leg  of  the  rectangle,  control  transfers  to  statement  50.  If  S is 
greater  than  0.95  times  BR2,  S is  set  equal  to  BR2.  The  x-  and  y-coordinates 

of  the  point  on  the  second  leg  are  computed  by  linear  interpolation  and  con- 

trol returns  to  the  calling  program. 

At  statement  50  the  x-  and  y-coordinates  of  a point  on  the  third 
* of  the  rectangle  are  computed  by  linear  interpolation.  Control  returns 
ie  calling  program. 

At  statement  60  the  distance,  S,  is  compared  to  the  length  of  the 
first  side  of  an  angle  segment.  If  S is  greater  than  this  length,  control 
transfers  to  statement  70.  If  not,  the  x-  and  y-coordinates  are  computed  by 
interpolation  for  a point  on  the  first  leg.  Control  returns  to  the  calling 
program. 


At  statement  70  the  distance  along  the  angle  is  decreased  by  the 
length  of  the  first  leg  of  the  angle.  The  coordinates  of  the  point  on  the 
second  leg  are  computed  by  linear  interpolation,  and  control  returns  to  the 
calling  program. 

i.  Subroutine  MAPPLT 

Subroutine  MAPPLT  draws  a map  of  the  street  segments,  one  line  per 
segment,  with  the  section  number  appended  to  each  segment.  Up  to  10  maps 
can  be  drawn. 

Subroutine  MAPPLT  has  two  arguments.  Argument  II  Indicates  the 
sequence  number  of  the  map.  Argument  KF  is  the  number  of  segments. 
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The  coordinates  of  the  region  bounding  the  map  are  contained  in 
arrays  in  COMMON  block  MPDATA.  In  this  COMMON  block,  arrays  XMIN  and  XMAX 
are  the  minimum  and  maximum  x-coordi nates  for  the  map.  XLEN  is  the  length, 
in  inches,  of  the  map  in  the  x-direction.  YMIN,  YMAX,  and  YLEN  are  the  corre- 
sponding arrays  in  the  y-direction.  Array  YHCUT  contains  the  height,  in  plot- 
ter inches,  at  which  the  map  must  be  sliced  into  strips.  Variable  AVMD  con- 
tains the  miles  per  MCU  conversion  factor  for  each  map. 

MAPPLT  begins  by  retrieving  or  computing  the  map  bounds,  the  height 
of  a strip  of  the  map  (PHGT),  the  maximum  length,  the  number  of  map  strips  (MX), 
the  map  scale  factors,  and  the  intervals  in  MCU  at  which  the  strips  are  to  be 
cut.  These  parameters  are  printed  according  to  format  90. 

The  loop  through  statement  200  will  test  each  segment  to  see  whether 
it  falls  within  the  frame  of  the  map;  if  it  does,  the  segment  will  be  plotted. 
Variables  NI  and  NF  are  set  equal  to  the  numbers  of  the  nodes  bounding  the 
segment.  The  midpoint  coordinates  of  the  segment  are  saved  in  variables  XMD 
and  YMD.  The  lines  in  the  node-number  array  at  which  the  initial  and  final 
nodes  occur  are  saved  in  variables  NS1  and  NS2.  The  initial  and  final  coor- 
dinates of  each  node  are  retrieved. 

Initially  the  segment  is  assumed  to  be  entirely  within  the  bounds, 
and  indicators  INBI,  INBM,  and  INBF  are  set  to  1 . If  the  coordinates  of  the 
initial  node  lie  outside  the  frame  of  the  map,  INBI  is  set  to  0.  Similar 
tests  are  made  on  the  coordinates  of  the  midpoint  of  the  segment  and  the  coor- 
dinates of  the  final  node  of  the  segment.  If  all  three  points  are  outside  the 
frame  of  the  map,  control  transfers  to  statement  200  and  the  segment  is  not 
plotted.  For  segments  that  are  at  least  partially  within  the  frame  of  the  map, 
the  section  number  and  total  length  of  the  segment,  in  miles,  are  saved  in 
variables  NUMS  and  T0TLEN.  The  number  of  points  to  be  used  in  plotting  one- 
half  of  the  segment  (NPMID)  is  computed.  The  number  will  be  restricted  to 
a maximum  of  10  points.  The  total  number  of  points  per  segment,  NPPSEG,  is 
set  to  twice  (NPMID). 

Subroutine  SHAPC0M  is  called  to  set  up  the  parameters  needed  to  gen- 
erate coordinates  of  points  on  the  segment.  The  cumulative  length  along  the 
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segment  is  initially  set  to  0.  A step  size,  DS,  is  computed  as  the  total 
length  divided  by  the  number  of  points  to  be  plotted  on  the  segment.  The  coor- 
dinates of  the  initial  node  are  stored  in  variables  XX  and  YY.  The  number  of 
the  strip  of  the  map  into  which  the  node  falls  is  computed.  Both  a current 
value  of  the  strip  number,  NMAP,  and  a value  for  the  previous  point,  NMAPO, 
will  be  used.  The  pen  position,  up  or  down,  is  determined  by  whether  the  ini- 
tial point  was  In  bounds.  Variable  IPEN  will  be  3 if  the  point  is  out  of  bounds 
and  2 if  the  point  is  in  bounds.  If  the  point  is  out  of  bounds,  control  trans- 
fers to  statement  130.  If  not,  the  coordinates  of  the  point  are  converted  to 
plotter  inches  and  stored  in  variables  XP  and  YP.  If  the  current  node  has  al- 
ready been  plotted  as  the  last  node  on  the  previous  segment,  control  transfers 
to  statement  120.  If  not,  a small  square  marking  its  position  is  appended  to 
the  map.  At  statement  120  the  pen  is  moved  to  the  position  of  the  current 
point  on  the  segment. 

Statement  130  starts  a loop  through  statement  170  that  will  advance 
the  pen  through  the  remaining  points  on  the  segment.  The  cumulative  length  is 
incremented  by  DS.  Subroutine  COORD  is  called  to  obtain  the  coordinates  of 
the  point  in  MCU. 

At  statement  140  the  coordinates  are  converted  to  plotter  inches. 

The  point  is  assumed  to  be  in  bounds,  and  variable  INB  is  set  to  1 . If  the 
coordinates  of  the  point  are  out  of  bounds,  INB  is  reset  to  0.  If  the  pen 
has  been  up  and  the  current  point  is  out  of  bounds,  or  if  the  strip  number 
is  greater  than  the  number  of  the  final  strip,  control  transfers  to  state- 
ment 60.  Otherwise,  the  pen  is  moved  to  the  position  of  the  current  point. 

If  the  pen  is  up,  it  is  lowered.  Variable  IPEN  is  recomputed  to  reflect 
whether  the  point  is  in  bounds. 

At  statement  150,  if  the  loop  index  is  not  equal  to  the  number  of 
the  midpoint  of  the  segment,  control  transfers  to  statement  160.  Otherwise, 
the  section  number  is  appended  to  the  map  near  the  segment  midpoint,  and  the 
pen  is  repositioned  at  the  midpoint. 

At  statement  160  the  number  of  the  current  strip  is  computed.  If 
the  current  strip  number  is  equal  to  the  previous  strip  number,  control  trans- 
fers to  statement  170.  If  not,  the  old  strip  number  (NMAPO)  is  set  equal  to 


22 


the  current  strip  number;  IPEN  is  set  to  3,  indicating  that  the  pen  is  up;  and 
control  transfers  to  statement  140.  In  this  case,  the  pen  is  positioned  at 
the  current  point  on  the  new  strip. 

Statement  170  is  the  end  of  the  loop  that  causes  the  segment  to  be 
drawn.  If  the  last  point  drawn  is  out  of  bounds,  control  transfers  to  state- 
ment 200.  Otherwise,  a small  square  marking  the  node's  position  is  appended 
to  the  map.  The  pen  is  repositioned  at  the  last  node.  The  number  of  the  node 
is  saved  in  variable  LASTNN.  Statement  200  is  the  end  of  the  loop  that  draws 
the  various  segments.  At  statement  300  the  plotter  pen  is  positioned  2 
inches  beyond  the  end  of  the  last  strip.  Control  returns  to  the  calling  pro- 
gram. 


j.  Subroutine  BUILD 

Subroutine  BUILD  creates  a near-neighbor  table  for  the  street  seg- 
ments in  the  map  description.  This  subroutine  has  13  arguments.  The  first,  N, 
is  the  total  number  of  segments  in  the  map  description.  The  second , KN,  is  the 
number  of  near  neighbors  that  will  be  found  for  each  segment.  KN  is  set  to  60 
in  the  main  program.  The  next  two  arguments,  X and  Y,  are  arrays  containing 
the  x-  and  y-coordi nates  of  the  segment  midpoints.  The  fifth  argument,  MINFR, 
is  an  array  containing  refuse  quantity,  servicing  time,  and  number  of  houses 
for  each  segment.  The  sixth  argument,  TREE,  is  an  array  used  in  the  construc- 
tion of  the  near-neighbor  table.  The  seventh  argument,  ISTPR,  is  an  array  of 
segment  numbers.  The  eighth  and  ninth  arguments,  NNT  and  NNTEMP,  are  tempor- 
ary storage  arrays.  The  next  two  arguments,  XT  and  YT,  are  arrays  of  the  x- 
and  y-coordi nates  of  the  segment  midpoints.  The  twelfth  argument,  KP,  is  the 
number  of  words  used  to  store  near-neighbor  information  for  each  segment.  The 
last  argument,  IUNX,  gives  the  number  of  the  unit  on  which  the  near-neighbor 
table  will  be  written.  If  IUNX  is  zero,  the  near-neighbor  table  will  be 
written  on  file  TAPE7. 

Subroutine  BUILD  begins  by  setting  variable  IUNIT  equal  to  UNIT(5). 
If  argument  IUNX  is  positive,  IUNIT  is  reset  to  IUNX.  The  loop  on  statement  5 
creates  sixty  1-bit  masks  in  array  BDATA.  The  next  statement  begins  a loop 
through  statement  1111  that  will  examine  each  segment.  The  segment  number  is 
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stored  in  NNT(l).  The  loop  through  statement  1000  transfers  the  segment  numbers 
to  array  NNTEMP  and  the  segment  midpoint  coordinates  to  arrays  XT  and  Y1 . The 
next  nine  statements  interchange  the  segment  stored  in  the  first  location  of 
arrays  XT,  YT,  and  NNTEMP  with  the  segment  stored  at  location  II.  Variable  MM 
is  set  to  one  less  than  the  number  of  segments.  The  loop  through  statement  20 
computes  the  distances  from  the  segment  whose  midpoint  coordinates  are  in  the 
first  location  of  arrays  XT  and  YT  to  each  other  segment.  The  distances  are 
stored  In  array  TREE  with  the  low-order  12  bits  replaced  by  the  segment  number. 
Subroutine  SORTK  Is  called  to  return  the  smallest  KN  distances  of  the  MM  dis- 
tances in  array  TREE.  The  loop  through  statement  30  retrieves  the  segment  num- 
bers from  the  low-order  12  bits  of  the  KN  smallest  distances.  These  segment 
numbers  are  stored  in  array  NNT.  The  loop  on  statement  94  sets  the  COMP  array 
equal  to  0.  This  array  will  be  used  to  store  in  turn  each  segment's  near- 
neighbor list.  The  loop  through  statement  95  generates  the  near-neighbor  list 
for  the  segment  currently  in  location  1 of  array  NNTEMP.  A unique  word  pointer, 
IW1,  and  a unique  bit  pointer,  I P 1 , are  generated  from  the  neighboring  segment 
number  in  array  NNT.  The  appropriate  bit  in  the  appropriate  word  of  array  COMP 
is  set  to  1 by  means  of  the  appropriate  mask  in  array  BDATA.  The  segment  num- 
ber is  stored  in  STRING(l).  The  load,  time,  and  number  of  houses  on  the  seg- 
ment are  transferred  to  array  MINFO,  which  is  equivalent  to  STRING(2).  Array 
STRING  is  written  to  file  IUNIT.  Note  that  the  COMP  array  is  equivalenced  to 
the  STRING  array  starting  at  STRING(5).  When  all  near-neighbor  table  informa- 
tion has  been  written  to  file  IUNIT,  the  file  is  rewound.  Control  returns  to 
the  calling  program. 

k.  Subroutine  SECTION 

Subroutine  SECTION  assigns  each  segment  in  the  map  description  to  a 
section  (a  section  corresponds  to  a collection  trip  or  vehicle  load).  Sub- 
routine SECTION  has  15  arguments.  The  first,  NN,  is  the  number  of  segments. 

The  second,  K,  is  the  number  of  near  neighbors  found  for  each  segment.  The 
next  two  arguments,  MODE  and  IFLAG,  are  provided  to  facilitate  modifications 
that  will  allow  for  user-specified  base  segments.  The  fifth  argument,  KCUT0F, 
gives  the  minimum  number  of  segments  to  be  considered  for  inclusion  in  the  same 
section  as  a given  base  segment.  The  sixth  and  seventh  arguments,  X and  Y,  are 
the  coordinates  of  the  segment  midpoints.  The  eighth,  NNTS,  is  an  array  of 
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segment  numbers.  The  next  four  arguments,  ISTO,  IST1,  1ST2,  and  IST4,  are  tem- 
porary storage  arrays  of  30  words  each.  The  thirteenth,  KP,  is  the  number  of 
words  used  for  near-neighbor  table  data  for  eudi  segment.  The  fourteenth,  KPB, 
is  the  number  of  words  in  use  per  segment  in  array  STRING.  The  fifteenth  argu- 
ment, MA,  is  the  dimension  of  arrays  ISTO,  IST1,  IST2,  and  IST4. 

Subroutine  SECTION  begins  by  initializing  parameters.  The  cumulative 
time  and  load  are  set  to  0.  The  maximum  number  of  base  segments  to  be  saved, 
NSTO,  is  set  to  30.  A count  of  the  number  of  passes  through  near-neighbor 
histogram  generation  (LPASS)  and  a count  of  the  number  of  completed  sections 
(KPASS)  are  set  to  0.  The  number  of  segments  is  stored  in  variable  NP.  Vari- 
able SMLD,  the  smal’est  cumulative  load  required  to  complete  the  current  sec- 
tion, is  set  to  0.  The  number  of  trucks,  NTRUCK,  is  set  to  0.  The  number  of 
nodes  is  saved  in  variable  N.  Symbolic  names  for  disk  files  are  assigned 
values.  Variable  0LDUNT  is  set  to  1,  IUNIT  to  2,  103  to  3,  104  to  4,  105  to 
7,  and  1010  to  10. 

The  loop  on  statement  5 clears  array  TRUCK,  while  the  loop  through 
statement  6 sets  the  cumulative  time  to  the  unloading  time,  TDUMP.  After 
statement  6,  the  pointer  to  the  first  segment  in  section  1 is  set  to  1.  The 
loop  through  statement  15  selects  vehicles  of  differing  capacities,  while  the 
loop  thorugh  statement  10  generates  an  entry  in  the  TRUCK  array  for  each  vehi- 
cle requested.  Following  statement  15,  the  total  number  of  vehicles  is  stored 
in  variable  NT0. 

Since  MODE  is  set  to  0 in  the  calling  program,  control  continues  to 
the  loop  through  statement  17.  (If  user-specified  base  segments  were  to  be 
added  to  the  program,  SECTION  would  be  called  with  a nonzero  MODE.)  The  loop 
through  statement  17  scans  the  near-neighbor  data  on  unit  7 (symbolically  105) 
searching  for  the  first  base  segment.  When  the  segment  Is  found,  the  segment 
data  and  the  neighbor  list  for  the  segment  are  transferred  to  array  BASE  in 
the  loop  on  statement  172.  The  refuse  quantity,  servicing  time,  and  number 
of  houses  for  the  segment  are  also  transferred  to  variables  INF1 , INF2,  and 
INF3.  Segments  other  than  the  base  segment  are  written  to  file  0L0UNT  at 
statement  171.  After  the  loop  through  statement  17  has  been  completed,  files 
0LDUNT  and  105  are  rewound.  Control  transfers  to  statement  1301,  bypassing 
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the  coding  that  selects  a base  segment  on  the  basis  of  its  distance  from  the 
previous  base  segment. 

At  statement  100  the  coordinates  of  the  current  base  segment  are 
transferred  to  variables  XR  and  YR.  A distance  variable,  ODIS,  is  set  to 
1000000.  The  loop  through  statement  1101  scans  the  segments  on  file  OLDUNT, 
computing  the  distance  from  each  segment  to  the  current  base  segment.  As  dis- 
tances shorter  than  ODIS  are  encountered,  the  new  distance  and  segment  numbers 
are  saved.  The  neighbor  data  are  saved  in  array  BASE.  When  the  loop  through 
statement  1101  is  complete,  the  segment  closest  to  the  old  base  segment  will 
be  saved  as  the  new  base  segment.  Files  OLDUNT  and  IUN IT  are  rewound. 

The  loop  through  statement  1401  scans  file  IUNIT  for  the  current 
base  segment.  All  other  segments  are  rewritten  to  file  OLDUNT.  When  the 
loop  through  statement  1401  is  complete,  both  files  are  rewound. 

Following  statement  1301,  the  smallest  cumulative  load  necessary  to 
complete  the  section  is  computed.  The  loop  on  statement  90  transfers  base  seg- 
ment information  to  array  NNTS.  Array  BASF,  contains  the  segment  number,  refuse 

quantity,  servicing  time,  number  of  houses,  and  near-neighbor  list.  This  infor- 
mation is  written  to  file  103.  The  count  of  the  number  of  base  segments  stored 
in  the  NNTS  array  (NSTD)  is  set  to  1.  The  base  segment  refuse  quantity,  servic- 
ing time,  and  number  of  houses  are  transferred  to  the  TRUCK  array.  Variable 
TRUCK(6,  SECTN)  is  set  to  1,  indicating  that  one  segment  is  serviced  by  the  ve- 
hicle in  this  section.  A count  of  segments  is  also  kept  in  variable  PC,  which 

is  also  set  to  1 . 

At  statement  1021  the  count  of  segments  left  to  be  assigned  (KL)  is 
computed.  At  statement  440,  counters  JON  and  ION  are  set  to  0.  The  loop 
through  statement  1020  sets  array  IST2  equal  to  consecutive  integers  and 
clears  arrays  IST0,  I ST  I , and  IST4.  The  loop  on  statement  1019  sets  to  0 the 
60  words  in  array  HIST0. 

The  segment  number  of  the  current  base  segment  is  stored  in  variable 
LI.  The  loop  through  statement  2929  scans  all  segments  on  file  OLDUNT.  The 
segment  and  its  neighbor  data  are  read,  and  the  segment  number  is  stored  in 
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variable  L2.  Word  pointer  IW1  and  bit  pointer  IP1,  which  indicate  the  position 
of  segment  L2  In  the  near-neighbor  data  for  the  base  segment,  are  computed.  This 
bit  Is  examined  in  the  neighbor  table  of  the  segment  read  from  OLDUNT;  If  it  Is 
nonzero,  control  transfers  to  statement  1022,  indicating  that  the  base  segment 
is  a near  neighbor  of  the  segment  from  file  OLDUNT.  Otherwise,  control  trans- 
fers to  statement  3030. 

At  statement  1022,  word  and  bit  pointers  are  computed  for  the  segment 
read  from  file  OLDUNT.  The  near-neighbor  table  for  the  base  segment  is  exam- 
ined. If  the  segment  from  file  OLDUNT  is  a near  neighbor  of  the  base  segment, 
control  transfers  to  statement  1023.  Otherwise,  control  transfers  to  statement 
3030. 

At  statement  1023,  the  shared-neighbor  count,  IC,  Is  set  to  0.  The 
loop  on  statement  1024  counts  the  number  of  neighbors  shared  by  the  base  seg- 
ment and  the  segment  read  from  file  OLDUNT.  The  count  is  used  as  a subscript 
on  array  HISTO,  and  the  appropriate  location  is  incremented  by  1 . If  the 
count  is  negative,  which  should  be  impossible,  or  if  variable  JON,  the  number 
of  segments  sharing  neighbors  with  the  base  segment,  is  greater  than  or  equal 
to  30,  control  transfers  to  statement  3030.  Otherwise,  JON  is  incremented  by 
1,  and  the  segment  number  and  shared-neighbor  count  are  saved  in  arrays  ISTO, 

I ST  1 , and  IST4.  The  STRING  array  for  this  segment  is  written  to  file  1010, 
and  control  transfers  to  the  end  of  the  loop. 

At  statement  3030,  the  STRING  array  for  segments  that  are  not  near 
neighbors  of  the  base  segment  is  written  to  file  IUNIT.  The  number  of  segments 
on  file  IUNIT  is  computed  and  stored  in  variable  NUT.  File  IUNIT  Is  rewound. 

The  loop  through  statement  4040  forms  a running  count  of  the  entries 
in  the  HISTO  array,  beginning  with  the  end  of  the  array.  When  the  count 
passes  KCUTOF,  control  transfers  to  statement  4450.  KCUTOF  is  set  to  5 in  the 
main  program.  At  least  five  segments  sharing  the  greatest  number  of  near 
neighbors  with  the  base  segment  will  be  examined  for  Inclusion  in  the  current 
section.  Variable  KI  is  set  equal  to  the  smallest  number  of  near  neighbors 
that  any  of  these  five  segments  shares  with  the  base  segment. 
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Following  statement  4450,  file  OLDUNT  is  rewound.  If  JON,  the  num- 
ber of  segments  sharing  neighbors  with  the  base  segment,  is  0,  control  trans- 
fers to  statement  4990.  Otherwise,  the  shared-neighbor  counts  in  array  IST1 
are  sorted  Into  decreasing  order;  the  segment  line  numbers  in  IST2  are  carried 
along  during  the  sort.  File  1010  is  rewound. 

A segment  counter,  variable  ION,  is  set  to  1 . At  statement  1001,  if 
ICODE  is  equal  to  1,  indicating  that  the  section  is  complete,  control  transfers 
to  statement  700.  Otherwise,  at  statement  101  variable  IP  is  set  equal  to  the 
line  number  of  the  unassigned  segment  sharing  the  most  neighbors  with  the  base 
segment.  Unit  1010  is  rewound.  The  segment  number  is  stored  in  variable  LNX, 
and  the  count  of  shared  neighbors  is  stored  in  LNY.  If  the  segment  counter 
(ION)  is  greater  than  the  total  number  of  segments  (JON),  control  transfers  to 
statement  4991.  Otherwise,  the  loop  on  statement  102  reads  segments  from  unit 
1010  until  a segment  is  found  with  a segment  number  equal  to  LNX  and  a shared- 
neighbor  count  greater  than  or  equal  to  variable  KI.  If  the  segment  is  found, 
control  transfers  to  statement  500.  If  not,  following  the  loop  through  state- 
ment 102,  if  the  current  load  plus  the  load  from  all  previous  sections  exceeds 
the  smallest  load  necessary  to  complete  the  current  section,  control  transfers 
to  statement  600.  Otherwise,  files  1010  and  OLDUNT  are  rewound. 

At  statement  4991,  counter  LPASS  is  incremented  by  1.  The  loop 
through  statement  3436  reads  the  segment  and  neighbor  data  from  file  1010.  If 
segments  have  not  been  assigned  to  the  current  section,  they  are  written  to 
file  OLDUNT.  The  loop  through  statement  3437  transfers  the  remaining  segments 
from  file  IUNIT  to  file  OLDUNT.  Files  IUNIT,  0LDL,. . , and  1010  are  rewound. 

At  statement  4990,  if  more  base  segments  are  needed  than  have  been 
saved  in  array  NNTS,  control  transfers  to  statement  800.  Otherwise,  the  loop 
on  statement  475  transfers  segment  and  neighbor  data  from  array  NNTS  to  array 
BASE.  Variable  NEXTN  is  incremented  by  1,  and  control  transfers  to  statement 
1021. 


At  statement  500  the  refuse  quantity  for  the  segment  is  stored  in 
variable  CURL.  The  servicing  time  on  the  segment  is  stored  in  variable  CURT. 
The  segment  counter,  ION,  is  incremented  by  1 . If  adding  the  segment  to  the 
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current  section  keeps  the  section  within  the  vehicle  capacity  and  time  limit, 
control  transfers  to  statement  550.  If  ION  is  greater  than  the  number  of  seg- 
ments sharing  neighbors  with  the  base  segment  (JON),  control  transfers  to  state- 
ment 600.  Otherwise,  file  1010  is  rewound.  Variable  IP  is  set  to  the  line  num- 
ber of  the  next  segment  sharing  neighbors  with  the  base  segment.  The  segment 
number  is  stored  in  variable  LNX,  and  the  count  of  shared  neighbors  is  stored 
in  variable  LNY. 

The  loop  through  statement  510  reads  the  segment  and  neighbor  data 
from  file  1010  into  array  STRING.  When  the  segment  being  sought  is  found,  con- 
trol transfers  to  statement  500.  If  the  segment  is  not  located  before  the  loop 
is  completed,  control  transfers  to  statement  600. 

At  statement  550  the  count  of  shared  neighbors  in  array  IST1  is  made 
negative,  indicating  that  the  segment  has  been  assigned  to  a section.  The  next 
three  statements  add  the  refuse  quantity,  the  servicing  time,  and  the  number  of 
houses  on  the  segment  to  the  corresponding  quantities  for  the  section  in  prog- 
ress. If  the  number  of  segments  saved  for  use  as  base  segments,  NSTD,  is 
greater  than  or  equal  to  the  maximum  number  allowable,  NST0,  control  transfers 
to  statement  598.  Otherwise,  NSTD  is  incremented  by  1. 

The  loop  on  statement  599  moves  the  current  segment  information  in 
array  STRING  to  the  base  segment  array  NSTS.  Following  statement  598,  the 
STRING  array  for  the  current  segment  is  written  to  file  103.  A segment  count, 
PC,  is  incremented  by  1.  The  number  of  segments  in  the  section  is  incremented 
by  1.  Control  transfers  to  statement  1001. 

Statement  600  is  reached  when  a section  is  full  or  when  no  other 
segments  can  be  added  to  the  section.  When  the  section  is  complete,  vari- 
able IC0DE  is  set  to  1.  File  1010  is  rewound. 

The  loop  through  statement  698  reads  the  segment  and  neighbor  data 
from  file  1010  into  array  STRING.  If  any  segment  has  not  been  assigned  to 
a section,  IST1  will  be  positive  for  that  segment  and  array  STRING  will  be 
written  to  file  0LDUNT.  After  the  loop  has  been  completed,  file  1010  is  re- 
wound. 
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The  loop  through  statement  699  transfers  the  remaining  segments  and 
neighbor  lists  from  file  IUNIT  to  file  OLDUNT.  Both  files  are  rewound  when 
the  loop  has  been  completed. 

At  stati'inent  700,  file  103  is  rewound.  The  count  of  total  unassigned 
segment  N,  is  decremented  by  PC,  the  count  of  segments  in  the  section  just 
completed.  Variable  NP  is  set  to  the  number  of  currently  unassigned  segments. 
Variable  LC  is  set  equa  to  PC.  The  loop  through  statement  705  reads  from  file 
103  the  segments  assigned  to  the  section  just  completed  and  adds  the  refuse 
quantity  and  servicing  time  to  the  cumulative  totals.  The  segment  numbers  are 
written  to  file  104. 

After  the  loop  through  statement  705  has  been  completed,  the  refuse 
quantity,  IESIL,  and  the  servicing  time,  TESTT,  for  all  remaining  unassigned 
segment  , are  computed.  The  number  of  the  next  section,  1ST,  is  computed.  If 
1ST  is  ireater  than  the  total  number  of  vehicles,  NTRUCK,  control  transfers 
to  statement  801. 

At  statement  818,  if  the  unass i lined  refuse  quantity  exceeds  the  ca- 
pacity >f  the  next  vehicle,  control  transfers  to  statement  710.  If  not,  and 
if  the  servicing  time  for  all  remaining  unassigned  segments  exceeds  the  time 
limit  for  the  next  vehicle,  control  transfers  to  statement  710.  Otherwise, 
the  count  of  total  segments  assigned,  PK,  is  incremented  by  PC.  The  sequence 
number  of  the  next  segment  to  be  assigned  is  stored  in  the  TRUCK  array.  File 
IUNIT  is  rewound.  Variable  LN  is  set  to  the  number  of  segments  yet  to  be 
assigned.  The  loop  through  statement  706  reads  the  remaining  unassigned  seg- 
ments from  file  OLDUNT  and  writes  the  segment  numbers  on  file  104.  When  the 
loop  is  complete,  the  section  count,  SECTN,  is  incremented.  Pointers  to  the 
first  and  last  segment  on  file  104  in  the  current  section  are  stored  in  array 
TRUCK.  File  104  is  rewound,  and  control  returns  to  the  calling  program. 

Statement  801  is  reacht  when  additional  trucks  must  be  defined  to 
complete  the  sectioning.  At  statement  801  the  number  of  vehicles,  NTRUCK,  is 
incremented.  A message  is  printed  indicating  that  the  vehicle  configuration 
has  been  extended.  Cyclical  counter  TPR  is  reset  to  1 if  it  exceeds  the  orig- 
inal number  of  vehicles.  A new  vehicle  is  defined  in  the  TRUCK  array  according 
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to  the  current  value  of  TPR.  The  loop  on  statement  8 11  clears  all  except  ca- 
pacity and  time-limit  items  in  the  TRUCK  array  on  the  line  for  the  new  vehicle. 
TPR  is  incremented  by  1.  Variable  1ST  is  set  equal  to  the  number  of  the  new 
vehicle.  Control  transfers  to  statement  818. 

Statement  710  is  reached  when  a section  is  completed.  At  statement 
710,  files  OLDUNT,  IUNIT,  and  103  are  rewound.  The  count  of  assigned  segments, 
PK,  is  incremented  by  PC.  The  section  number  is  incremented  by  1 . A pointer 
to  the  next  segment  to  be  written  to  file  104  is  computed  and  stored  in  the 
TRUCK  array.  Various  parameters  are  reset,  and  control  transfers  to  statement 
100. 


Statement  800  is  reached  when  base  segments  can  no  longer  be  obtained 
from  the  NNTS  array.  The  coordinates  of  the  last  base  segment  midpoint  are 
stored  in  variables  XR  and  YR.  A distance  variable,  0DIS,  is  initially  set  to 
1000000.  The  number  of  remaining  unassigned  segments  is  computed  and  stored 
in  variable  NLK.  The  count  of  base  segments  is  reset  to  0. 

The  loop  through  statement  2101  searches  for  a new  base  segment. 
Segments  are  read  from  file  OLDUNT  and  are  written  to  file  IUNIT.  The  segment 
number  is  saved  in  variable  NS.  The  segment  midpoint  coordinates  are  stored 
in  variables  XS  and  YS.  The  distance  between  the  current  segment  and  the  pre- 
vious base  segment  is  computed  and  stored  in  variable  DIS.  If  DIS  is  greater 
than  or  equal  to  0DIS,  control  transfers  to  the  end  of  the  loop.  If  the  refuse 
quantity  for  the  current  segment  causes  the  vehicle  capacity  to  be  exceeded, 
control  transfers  to  statement  2101.  If  not,  and  if  the  servicing  time  for 
the  current  segment  causes  the  vehicle  capacity  to  be  exceeded,  control  trans- 
fers to  statement  2101.  Otherwise,  the  current  distance  is  saved  in  0DIS,  and 
the  current  segment  number  is  saved  in  variable  NBASE. 

The  loop  on  statement  2101  transfers  the  STRING  array  for  the  seg- 
ment to  the  BASE  array.  Statement  2101  ends  the  loop  that  seeks  a new  base 
segment.  Files  OLD  and  IUNIT  are  rewound.  If  no  base  segment  was  found, 

NBASE  will  be  equal  to  0,  and  control  will  transfer  to  statement  600. 
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The  loop  through  statement  2401  reads  segments  and  their  neighbor 
lists  from  file  I UN  I T into  array  STRING.  If  the  segment  number  is  not  equal 
to  the  base  segment  number,  control  transfers  to  statement  2400.  Otherwise, 
the  segment  refuse  quantity,  servicing  time,  and  number  of  houses  are  trans- 
ferred to  variables  INF  1 , I NF2 , and  INF3.  Control  transfers  to  statement  2401. 

At  statement  2400  the  STRING  array  is  written  to  file  0LDUNT.  Fol- 
lowing the  loop,  files  0LDUNT  and  IUNIT  are  rewound.  The  base  segment  and  its 
neighbor  data  are  written  to  file  103.  The  count  of  segments  assigned  to  the 
current  section,  PC,  is  incremented  by  1.  Section  data  in  array  TRUCK  are  up- 
dated for  the  refuse  quantity,  servicing  time,  number  of  segments,  and  number 
of  houses.  The  number  of  base  segments,  NSTD,  is  reset  to  1 and  variable  NFXTN 
is  set  to  2.  If  the  total  amount  of  refuse  assigned  so  far  is  less  than  the 
smallest  load  required  to  complete  the  current  section,  control  transfers  to 
statement  1021.  Otherwise,  control  transfers  to  statement  600. 

1.  Program  PHASE2 

Main  program  PHASE2  drives  the  sectioning  and  the  plotting  of  section 
assignment  maps.  It  uses  11  files:  INPUT,  OUTPUT,  TAPE1 , TAPE2 , TAPE3 , TAPE4, 
TAPE7,  TAPES,  TAPE9,  TAPE10,  and  TAPE! 1 . File  TAPE5  is  equivalenced  to  INPUT 
in  order  to  test  for  end-of-record  cards  in  the  card  input  data.  File  TAPEl 
is  used  for  section-description  information.  Files  TAPE2  and  TAPE3  are  used 
for  temporary  storage  of  segment  and  near-neighbor  information.  File  TAPE4  is 
used  for  a list  of  segment  numbers  in  the  order  they  are  assigned  to  sections. 
File  TAPE7  is  used  for  temporary  storage  of  segment  and  neighbor  data.  File 
TAPE8  is  the  Calcomp  plot  tape.  File  TAPE9  holds  input  segment  data  from  pro- 
gram RCINPT.  File  TAPE10  is  used  for  temporary  storage  of  segment  and  neighbor 
data.  File  TAPE11  holds  node  data  from  program  RCINPT. 

Blank  COMMON  holds  arrays  for  the  problem  title  and  the  segment  data. 
Array  ISEG  contains  the  section  numbers  assigned  to  segments  by  subroutine 
SECTION.  Arrays  NN1  and  NN2  are  the  starting  and  ending  node  numbers  for  the 
segments.  Array  FLEN  holds  the  lengths  of  the  segments.  Array  NH  is  the  num- 
ber of  houses  on  the  segments.  Array  FMPH  is  the  speed  limits  on  the  segments. 
Array  RQF  gives  the  refuse  quantity  adjustment  factor  for  the  segments.  Arrays 
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X and  Y are  the  midpoint  coordinates  of  the  segments.  Array  SF  holds  the  shape 
codes  for  the  segments. 

COMMON  block  MPDATA  holds  arrays  describing  the  map  bounds  and  sizes 
for  up  to  ten  maps.  Array  YHCUT  contains  the  height  of  a strip  of  the  map,  and 
variable  AVMO  is  the  number  of  miles  per  MCU. 

COMMON  block  NDDATA  bolds  the  node  data.  KNODES  is  the  number  of 
nodes.  Array  NBS  holds  the  segment  numbers  of  up  to  six  segments  bounding 
each  node.  Array  NODNUM  contains  the  node  numbers.  Arrays  XNOD  and  YNOD  are 
the  coordinates  of  the  nodes. 

COMMON  block  DISKIO  holds  array  UNOT,  which  is  used  for  symbolic 
references  to  disk  and  tape  file  numbers. 

COMMON  block  ROUTE  holds  vehicle  and  section-description  data.  Array 
TRUCK  has  seven  words  of  data  for  each  of  up  to  50  vehicles.  Array  TRUCKS  holds 
descriptive  information  for  the  vehicle  fleet  available  in  the  program.  Vari- 
able NTRUCK  is  the  total  number  of  vehicles  or  sections.  Variable  NBASE  is  the 
segment  number  for  the  segment  to  be  used  as  the  first  base  segment. 

COMMON  block  STATS  holds  additional  information  about  the  sections. 
Variable  FRACT  is  the  ratio  of  total  refuse  to  total  vehicle  capacity.  Vari- 
ables TOTL  and  TOTT  are  the  total  refuse  quantity  and  total  servicing  time  for 
all  segments  in  the  map  description.  Variables  CUML  and  CUMT  are  the  cumula- 
tive load  and  servicing  time.  Variable  ISECTN  is  set  to  the  total  number  of 
sections  by  subroutine  SECTION.  Variable  TDUMP  is  the  unloading  time  at  the 
landfill . 


Program  PHASE2  begins  execution  by  reading  three  data  cards:  the 
title  card,  the  vehicle-description  card,  and  the  time-limits  card.  Segment 
data  are  read  from  file  TAPE9.  Node  data  are  read  from  file  TAPE11. 

The  loop  on  statement  15  retrieves  the  absolute  value  of  the  number 
of  houses  on  each  segment.  A negative  number  for  variable  NH  indicates  that 
collection  is  from  the  right  side  of  the  street  only.  The  number  of  segments 
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is  stored  in  variable  NA.  If  the  number  of  the  first  base  segemnt,  NBASE,  was 
left  blank  on  the  third  data  card  or  was  improperly  specified,  it  is  set  to  1. 
Variable  KN,  the  number  of  near  neighbors  to  be  found  for  each  segment,  is  set 
to  60.  Variables  MODE,  OPTION,  and  NSEEO,  which  are  provided  to  facilitate  im- 
plementation of  user-specified  base  segments,  are  set  to  0.  The  total  time  and 
load  are  set  to  0. 

The  loop  through  statement  20  calculates  and  saves  the  total  refuse 
quantity,  servicing  time,  and  number  of  houses  for  each  segment.  On  each  seg- 
ment with  houses,  the  servicing  time  is  recomputed  to  include  travel  at  5 miles 
per  hour  and  stopping  time  at  the  houses.  The  total  time  and  load  in  the  net- 
work is  accumulated  in  variables  TOTT  and  TOTL.  Following  statement  20,  if  the 
maximum  trip  time  is  unspecified  or  negative,  it  is  reset  to  24  hours.  The  fol- 
lowing statement  converts  the  maximum  trip  time  to  minutes.  For  a vehicle  with 
a nonzero  capacity,  the  loop  on  statement  25  sets  the  number  of  vehicles  in  the 
NT  array  to  1 if  no  vehicles  were  specified.  Variable  REF  is  initially  set  to 
0.  It  will  be  used  to  accumulate  the  total  capacity  available. 

The  loop  through  statement  50  moves  to  the  TRUCKS  array  the  number 
of  vehicles,  the  capacity,  and  the  maximum  trip  time  for  each  of  the  four  types 
of  vehicles  allowed.  The  total  vehicle  capacity  available,  REF,  is  also  accumu- 
lated in  this  loop.  The  integer  part  of  the  ratio  of  total  refuse  to  total  ve- 
hicle capacity  is  computed  and  stored  in  variable  NTEA.  If  NTEA  is  0,  which 
indicates  that  enough  vehicles  were  specified  in  the  card  input  to  completely 
service  the  region,  control  transfers  to  statement  80.  Otherwise,  REF  is  re- 
set to  0. 


The  loop  through  statement  60  multiplies  the  number  of  vehicles  spe- 
cified on  inDut  cards  by  NTEA  and  stores  the  result  in  the  TRUCKS  array.  The 
vehicle  capacity  now  available  is  accumulated  at  statement  60.  The  loop 
through  statement  70  adds  additional  vehicles  where  necessary  to  bring  the 
total  fleet  capacity  above  the  total  amount  of  refuse  to  be  collected.  An  in- 
crement, NINC,  is  preset  to  0.  For  each  vehicle  capacity  specified,  NINC  is 
set  to  the  smaller  of  the  number  of  vehicles  specified  on  the  data  card  or 
the  number  of  vehicles  required  to  provide  enough  capacity  to  completely  ser- 
vice the  region.  The  vehicle  capacity  available,  REF,  is  increased  by  NINC 
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times  the  vehicle  capacity.  The  number  of  vehicles  is  increased  by  NINC  in 
the  TRUCKS  array.  If  enough  capacity  is  available  to  service  the  region,  con- 
trol transfers  to  statement  80. 


At  statement  80  the  ratio  of  total  refuse  to  total  vehicle  capacity 
is  computed  and  stored  in  variable  FRACT.  The  input  card  data  are  printed, 
along  with  FRACT,  the  minimum  fraction  that  each  vehicle  must  be  filled.  The 
starting  base  segment,  NBASE , is  also  printed.  Variable  IUNX  is  set  to  0. 

The  loop  on  statement  96  generates  an  array  of  segment  numbers 
(ISTPR).  The  segments  are  numbered  sequentially  starting  at  1.  Variable  MA, 
the  maximum  number  of  base  segments  to  be  saved  while  building  a section,  is 
set  to  30.  Variable  KP,  the  number  of  words  needed  to  store  the  near-neighbor 
data,  is  computed.  Variable  KPB,  the  number  of  words  used  in  array  STRING,  is 
set  to  KP  plus  4.  Subroutine  BUILD  is  called  to  generate  the  near-neighbor 
table.  Variable  IFLAG  is  set  to  0.  Variable  KCUT0F  is  set  to  5,  causing  at 
least  five  segments  sharing  the  most  neighbors  with  a base  segment  to  be  ex- 
amined  for  addition  to  a section. 

Subroutine  SECTION  is  called  to  assign  the  segments  to  sections.  If 
IFLAG  is  set  to  1,  control  transfers  to  statement  333.  Otherwise,  a sumnary  is 
printed  for  the  various  trips.  The  summary  includes  trip  number,  vehicle  ca- 
pacity, vehicle  time  limit,  vehicle  load,  servicing  time  for  the  section,  num- 
ber of  segments  in  the  section,  and  number  of  houses  in  the  section.  The  loop 
through  statement  40  prints  the  data  in  this  tabulation.  If  OPTION  is  equal  to 
1,  the  number  of  segments  is  incremented  by  the  number  of  user-specif ied  base 
segments.  (Note  that  user-specified  base  segments  are  not  implemented  at  this 
time.)  Variable  IUNIT  is  set  to  UN0T(4),  which  has  value  4.  Variable  J is  set 
to  1. 

The  loop  through  statement  33  reads  the  segments  from  file  IUNIT  and 
groups  and  prints  them  according  to  their  section.  When  the  loop  index  is  not 
equal  to  the  sequence  number  of  a segment  starting  a section,  control  transfers 
to  statement  88.  Otherwise,  a heading  for  the  section  is  printed.  The  section 
number  J is  incremented  by  1 . A carriage  control,  variable  CC,  Is  set  to  a 
blank.  The  number  of  segments  on  the  current  line  of  printed  output,  NN,  is 
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set  to  0.  At  statement  88,  NN  is  incremented  by  1.  The  next  statement  prints 
the  segment  number  preceded  by  the  number  of  blanks  required  to  put  the  segment 
number  in  its  appropriate  position  on  the  line.  The  carriage  control  is  set  to 
a plus  sign.  If  fewer  than  30  segment  numbers  have  been  written,  control  trans- 
fers to  statement  33.  Otherwise,  the  carriage  control  is  reset  to  a blank,  and 
the  count  of  segment  numbers  on  the  line  is  reset  to  0.  Statement  33  ends  the 
loop  that  prints  the  segments  in  each  section.  File  IUNIT  is  rewound. 

The  plot  package  is  initialized  by  a call  to  PLOTS.  Subroutine  PLOT 
is  called  twice,  moving  the  pen  down  and  then  up  3 inches  in  order  to  create 
a 3-1nch  border  on  the  plot. 

The  loop  through  statement  1000  controls  the  reading  of  segments  by 
section.  Variable  IL  is  set  to  the  sequence  number  of  the  last  segment  in  the 
Ith  section.  The  loop  through  statement  2000  reads  each  segment  from  file 
IUNIT  and  assigns  its  section  number  to  the  appropriate  location  in  array  ISEG. 
The  number  of  maps,  MAPS,  is  initially  set  to  0.  The  loop  through  statement 
1030  controls  the  reading  of  up  to  10  output-map-bounds  cards.  The  minimum  and 
maximum  coordinates  and  lengths  in  the  x and  y direction  are  read  from  a card. 
Also  read  is  the  height  at  which  the  map  should  be  cut  into  strips.  If  an  end- 
of-file  card  is  encountered  during  the  read,  control  transfers  to  statement 
1040.  Otherwise,  execution  continues  at  statement  1020.  At  statement  1020, 
if  the  map-strip  height  is  less  than  or  equal  to  zero,  or  is  unspecified  on  the 
card,  the  strip  height  will  be  set  to  30  inches.  At  statement  1030  the  number 
of  maps  is  set  equal  to  the  loop  index,  I.  At  statement  1040,  if  the  number  of 
maps  is  greater  than  zero--that  is,  if  any  map-bojnds  cards  have  been  found-- 
control  transfers  to  1070.  Otherwise,  default  values  are  set. 

The  default  values  are  set  in  the  following  manner.  MAPS  is  set  to 
1,  and  the  lengths  in  the  x and  y directions  are  set  to  30  inches.  The  height 
of  a map  strip  is  set  to  30  inches.  Bounds  on  the  coordinates  are  set  ini- 
tially so  that  the  minimums  are  large  positive  numbers  and  the  maximums  are 
large  negative  numbers.  The  loop  through  statement  1050  scans  the  midpoint 
coordinates  of  the  segments  and  resets  the  minimum  and  maximum  coordinates 
appropriately.  The  loop  through  statement  1060  scans  the  node  coordinates 
and  resets  the  minimum  and  maximum  coordinates  appropriately. 
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The  loop  on  statement  1080  calls  subroutine  MAPPLT  to  plot  each  out- 
put map.  Subroutine  PLOT  is  called  to  terminate  the  plot  file.  Following 
statement  333,  file  TAPE1  is  rewound.  The  number  of  segments  and  the  number 
of  sections  are  written  to  TAPE1 , followed  by  the  sequence  numbers  of  the  first 
and  last  segments  in  each  section  and  the  vehicle  capacity.  Files  TAPE1  and 
TAPE4  are  end-filed,  and  program  execution  stops. 


(The  reverse  of  this  page  is  blank.) 
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SECTION  IV 
INPUT  AND  OUTPUT 

1.  INPUT 

Input  to  program  PHASE2  consists  of  card  input,  segment  data  from  file 
TAPE9,  and  node  data  from  file  TAPE  1 1 . 

a.  Card  Input 

Table  1 presents  the  form  and  contents  of  the  four  types  of  data 
cards.  The  first  card  contains  a title,  which  is  printed  on  the  first  line  of 
the  output.  The  second  card  contains  the  number  and  capacity  of  up  to  four 
kinds  of  vehicles.  The  third  card  contains  time  restrictions  and  the  number  of 
the  segment  to  be  used  as  a base  segment  for  the  first  section.  The  fourth  card 
specifies  coordinate  bounds  and  sizes  for  the  output  map.  If  the  output  map  Is 
to  be  plotted  in  strips  less  than  30  inches  high,  the  height  of  the  strip  must 
be  indicated  on  the  fourth  card.  The  fourth  card  may  be  omitted,  or  it  may  be 
repeated  up  to  10  times.  If  it  is  omitted,  one  30-  by  30-1nch  map  of  the  en- 
tire collection  region  will  be  plotted.  If  It  Is  repeated,  any  map-bounds  cards 
after  the  tenth  card  will  be  ignored. 

b.  Segment  Data 

Disk  file  TAPE9  contains  segment  data  and  the  map  distance  conversion 
factor  (miles  per  MCU)  for  the  overall  map.  All  of  the  data  are  read  by  one 
binary  READ  statement.  The  first  word  is  the  count  of  the  segments.  The  seg- 
ment data  follow,  11  words  per  segment  for  each  segment.  After  the  segment 
data  comes  the  overall  distance  conversion  factor. 

The  list  used  in  the  READ  statement  Is  NSEG,(DUMMY,NN1(I),NN2(I), 

FLEN ( I ) , NH ( I ) , FMPH ( I ) ,DUMMY,RQF( I ) ,X( I ) ,Y( I ) ,SF( I ) , 1=1 ,NSEG) ,AVMD.  In  the 
list,  variable  NSEG  is  the  count  of  segments.  Since  the  street  numbers  of  the 
segments  are  not  needed,  they  are  read  into  variable  DUMMY.  Arrays  NN1  and  NN2 
hold  the  starting  and  ending  node  numbers  for  the  segments.  Array  FLEN  holds 


39 


TABLE  1.  PHASE 2 DATA  CARDS 


61-70 


F10.0 


Height  of  map  strips.  In  Inches 


I 


I 
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the  length  of  the  segment,  in  miles.  Array  NH  holds  the  number  of  houses. 

Array  FMPH  holds  the  speed  limits  for  the  segments.  The  one-way-street 
indicators  are  not  needed  and  are  read  into  variable  DUMMY.  The  refuse  quan- 
tity adjustment  factors  are  read  into  array  RQF.  The  midpoint  coordinates  of 
the  segments  are  read  into  arrays  X and  Y.  The  shape  codes  for  the  segments 
are  read  into  array  SF.  Variable  AVMD  is  the  map-distance  conversion  factor. 

c.  Node  data 

Disk  file  TAPE11  contains  refuse-quantity  information  and  node  data. 

All  of  the  data  are  read  by  one  binary  READ  statement.  The  list  used  in  the 
READ  statement  is  NHTOT.TOTREF, KNODES, (NODNUM(I),NBS(I),XNOD(I),YNOD(I), 1=1, 

KNODES).  The  first  three  words,  NHTOT,  TOTREF,  and  KNODES,  are  the  total  num- 
ber of  houses,  the  total  refuse  quantity,  and  a count  of  the  nodes.  Array 
NODNUM  holds  the  node  numbers.  Array  NBS  holds  up  to  six  segment  numbers  of 
segments  bounding  the  node.  Arrays  XNOD  and  YNOD  are  the  coordinates  of  the 
nodes.  f 

2.  SCRATCH  FILES 

Disk  files  TAPE1,  TAPE2,  TAPE3,  TAPE7,  and  TAPE10  are  used  by  subroutine 
SECTION  as  scratch  files,  and  all  contain  the  same  type  of  data.  The  data  con- 
sist of  40  words  of  segment  and  near-neighbor  information  for  each  segment.  The 
first  word  is  the  segment  number.  The  second  word  is  the  refuse  quantity.  The 
third  word  is  the  servicing  time  for  the  segment.  The  fourth  word  is  the  number 
of  houses  on  the  segment.  The  next  25  words  are  the  near-neighbor  list  for  the 
segment.  The  last  11  words  are  not  used  at  present.  The  first  word  of  the 
near-neighbor  list  indicates  which  of  the  first  60  segments  are  near  neighbors 
to  the  segment.  The  low-order  bit  corresponds  to  segment  number  1.  The  bits 
are  set  to  1 If  the  corresponding  segment  is  one  of  the  60  nearest  segments; 
otherwise,  they  remain  0.  The  second  word  of  the  neighbor  list  contains  the 
bits  corresponding  to  segments  61  through  120,  and  so  forth. 

The  segment  data  and  the  neighbor  list  are  written  Initially  to  TAPE7  by 
subroutine  BUILD.  TAPE7  is  read  in  subroutine  SECTION,  and  all  segment  data 
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except  those  for  the  first  base  segment  are  rewritten  to  file  TAPE1,  symbol- 
ically referred  to  as  OLDUNT.  The  first  base  segment  is  written  to  file  TAPE3. 
TAPE7  is  not  used  again.  The  data  on  OLDUNT  are  read,  and  segments  sharing 
near  neighbors  with  the  base  segment  are  written  to  file  TAPE10,  symbolically 
1010.  Those  segments  not  sharing  neighbors  with  the  base  segment  are  written 
to  file  TAPE2,  symbolically  IUNIT.  The  data  on  TAPE10  are  examined,  and  those 
which  can  be  added  to  the  current  section  are  written  on  TAPE3.  Those  which 
cannot  be  added  are  written  to  OLDUNT.  The  segments  on  IUNIT  are  copied  to 
OLDUNT,  and  the  process  repeats  until  a section  is  completed.  When  a section 
is  completed,  the  segment  numbers  are  transferred  from  TAPE3  to  TAPE4.  When 
only  one  section  remains  to  be  completed,  all  segment  numbers  from  OLDUNT  are 
copied  directly  to  TAPE4. 

In  sunnary,  TAPE1  holds  the  segment  and  neighbor  data  for  segments  not  yet 
assigned  to  a section.  TAPE2  holds  segment  and  neighbor  data  for  all  segments 
not  sharing  enough  neighbors  with  the  base  segment.  TAPE3  holds  segment  and 
neighbor  data  for  segments  added  to  the  current  section.  TAPE7  holds  segment 
and  neighbor  data  for  all  of  the  segments.  TAPE10  holds  segment  and  neighbor 
data  for  those  segments  sharing  enough  neighbors  with  the  base  segment  to  be 
considered  for  addition  to  the  section. 


3.  OUTPUT 

a.  Disk  and  Plot  Files 

File  TAPE1  Is  reused  for  output  in  main  program  PHASE2  after  the  map 
plotting  Is  completed.  The  list  used  in  the  binary  WRITE  statement  Is  NA, 
ISECTN,(TRUCK(5,I),TRUCK(6,I)+TRUCK(5,I)-1 TRUCK( 1 ,1 ) ,1=1 .ISECTN) . Variable 
NA  is  the  number  of  segments;  ISECTN  is  the  number  of  sections.  The  next  three 
items  written  for  each  section  are  the  sequence  number  on  TAPE4  of  the  first 
segment  in  the  section,  ihe  sequence  number  of  the  last  segment  in  the  section, 
and  the  vehicle  capacity.  The  file  should  be  cataloged  for  later  use  by  pro- 
gram PHASE3. 
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Disk  file  TAPE4  is  used  to  hold  the  segment  numbers  In  the  order  in 

which  they  are  assigned  to  sections.  The  segment  numbers  are  written  with  a 

formatted  WRITE  statement,  one  segment  number  at  a time.  The  format  used  is 
IX, 15.  The  file  should  be  cataloged  for  later  use  by  program  PHASE3. 

File  TAPE8,  the  plot  file,  will  be  disk  or  tape  depending  on  the  pro- 
cedure used  by  the  local  installation  to  produce  plots.  Each  output  map  re- 
quested occupies  one  file.  PHASE2  generates  an  empty  file  before  terminating 
TAPE8.  The  structure  of  the  file  depends  on  the  local  installation.  Figure  2 
shows  a section  map  for  Kirtland  Air  Force  Base. 

b.  Printed  Output 

There  are  four  parts  to  the  printed  output:  a list  of  the  input  card 

data,  a table  of  trip  information,  a list  of  segments  In  each  section,  and  a 

list  of  the  parameters  used  in  each  output  map. 

Appendix  D contains  sample  printed  output.  The  first  page  summarizes 
the  vehicle  input  data.  The  title  is  printed  following  the  heading  INPUT 
VEHICLE  DATA.  The  capacity  and  number  of  vehicles  available  for  four  kinds  of 
vehicles  are  listed.  Zeros  are  printed  where  no  vehicle  was  specified.  If 
the  number  of  vehicles  on  the  input  data  was  inadequate  to  service  the  region, 
the  number  of  trips  required  will  be  listed  in  the  column  headed  TRIPS.  The 
ratio  of  total  refuse  to  total  vehicle  capacity  is  printed  as  MINIMUM  FILL 
FRACTION.  The  sectioning  algorithm  endeavors  to  fill  each  vehicle  to  the  in- 
dicated fraction  before  completing  a section.  The  next  four  lines  give  the 
times  specified  on  the  third  data  card.  The  segment  to  be  used  as  the  first 
base  segment  is  also  printed.  If  this  itez<  has  been  left  blank  on  the  third 
data  card,  'the  printed  message  indicates  that  sectioning  starts  with  segment  1. 

The  second  page  of  printed  output  starts  with  a table  of  Information 
about  each  section.  The  capacity,  time  limit,  and  load  for  the  vehicle  is 
given,  followed  by  the  collection-time  estimate,  the  number  of  segments,  and 
the  number  of  houses  in  the  section.  At  this  stage,  the  collection  time  is 
approximate  because  some  streets  not  requiring  collection  may  be  dropped  from 
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inciuuea  in  tne  co  iiecuon-iime  estimate. 


The  printed  output  continues  with  a list  of  the  segments  in  each 
section.  The  segment  numbers  are  listed  in  the  order  in  which  they  were  se- 
lected for  inclusion  in  the  section  and  will  not  usually  be  in  numerical 
order.  UP  to  30  segment  numbers  are  printed  on  each  line. 

Fourteen  parameters  are  printed  for  each  output  map,  following  the 
heading  "MAPPLT  PARAMETERS  FOR  MAP  n,"  where  n is  the  sequence  number  of  the 
map.  AVMD  is  the  map  distance  conversion  factor  in  miles  per  MCU.  Each  of 
the  next  two  lines  gives  the  minimum  and  maximum  x-  and  y-coordinates  In  MCU. 

On  the  final  line,  XSC  and  YSC  identify  the  x and  y map  scales  in  Inches  per 
MCU.  The  plot  strip  height,  PHGT,  and  the  overall  plot  length,  PLEN,  are 
printed  in  inches.  YCUT  is  the  height  of  a map  strip  In  MCU. 

V 
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SECTION  V 

PROGRAM  REQUIREMENTS 


1 . SYSTEM 

Except  for  function  KOUNT,  program  PHASE2  is  written  entirely  in  FORTRAN 
IV.  Function  KOUNT  is  written  in  the  COMPASS  assembler  language.  The  program 
runs  on  a CDC  6600  using  a SCOPE  3.4.4  operating  system. 

Eleven  obvious  types  of  computer-dependent  coding  are  used  in  program 
PHASE2  and  its  subroutines.  A 60-bit  word  is  assumed  in  the  main  program  and 
In  subroutines  BUILD,  SECTION,  and  NUMBER.  System  function  SHIFT  is  used  In 
subroutine  BUILD.  An  R format  is  used  in  subroutines  SHAPCOM  and  COORD.  Six- 
bit  characters  are  assumed  in  subroutine  NUMBER.  An  ENCODE  statement  Is  used 
in  subroutine  NUMBER.  In-line  function  MINO  Is  used  in  PHASE2.  Asterisk- 
bounded  character  strings  are  used  in  formats  in  PHASE2  and  in  subroutine 
SECTION.  A computation  is  used  as  an  item  in  an  output  list  in  PHASE2.  The 
AND  operation  is  used  for  masking  in  subroutine  SECTION.  Function  KOUNT  is 
written  in  COMPASS.  A dollar  sign  is  used  between  FORTRAN  statements  in  pro- 
gram PHASE2  and  subroutines  SORTK,  SECTION,  MAPPLT,  SHAPCOM,  and  COORD.  More 
subtle  types  of  machine  dependencies  may  exist,  according  to  the  machine  used. 


2.  STORAGE 

The  core  requirement  is  slightly  less  than  114,0008  words,  but  may  vary 
slightly  depending  on  the  plotting  system  used  by  the  local  Installation.  The 
peripheral  storage  for  output  disk  files  should  not  exceed  28,000  words  for 
TAPE1 , TAPE2,  TAPE3,  TAPE7,  and  TAPE10.  When  file  TAPE1  Is  cataloged.  It  should 
not  contain  more  than  152  words.  Disk  file  TAPE4  should  not  exceed  700  words. 
The  plot  file,  TAPE8,  should  not  exceed  1.5  million  words,  although  more  typ- 
ically the  file  will  contain  roughly  20,000  words  per  output  map. 
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3.  TIME 


The  execution  CP  time  for  program  PHASE2  Is  approximately  0.005  SjgT  + 
0.004  Jjps  S.  seconds,  where  STQT  Is  the  total  number  of  segments  and  is  the 
number  of  segments  on  the  1^  output  map.  The  maximum  possible  CP  time  using 
10  full  maps  of  a 700-segment  network  would  be  273  seconds.  Rough  estimates  of 
the  10  and  PP  times  are  0.2  Sjqt  seconds  for  10  time  and  0.5  STqT  seconds  for 


PP  time. 
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SECTION  VI 
PROGRAM  LIMITATIONS 


1 .  PLOTTER 

The  plotter  used  by  PHASE2  can  be  a drum  plotter  with  at  least  a 10-inch 
drum  or  a flatbed  plotter  with  a 30- inch-square  or  larger  bed.  If  the  user 
does  not  indicate  a plot-strip  height  on  the  fourth  data  card,  the  program  as- 
sumes that  a 30-inch  height  is  available.  Maps  may  not  be  drawn  in  strips  on 
a flatbed  plotter;  therefore,  the  height  and  width  of  the  output  map  must  not 
be  specified  larger  than  the  plot-strip  height.  Some  systems  limit  plot 
lengths;  for  example,  Calcomp  plots  may  not  exceed  120  Inches  in  length  at  the 
Air  Force  Weapons  Laboratory  at  Kirtland  Air  Force  Base.  Each  strip  of  an  out- 
put map  has  a 1-inch  space  after  it,  with  an  additional  inch  after  the  last 
strip.  Thus  a 30-  by  30-inch  map  plotted  as  three  30-  by  10-Inch  strips  gen- 
erates a plot  94  inches  long.  Any  limitations  on  plot  lengths  imposed  by  the 
local  system  must  be  considered  when  the  output-map  description  cards  are 
punched. 


2.  NODE  AND  SEGMENT 

Array  dimensions  in  program  PHASE2  limit  the  number  of  nodes  to  500  and 
the  number  of  segments  to  700.  Since  these  data  are  passed  to  PHASE2  by  pro- 
gram RCINPT  as  disk  files  TAPE9  and  TAPE11,  the  limits  should  not  be  exceeded. 


3.  VEHICLE  FLEET 

The  vehicle  fleet  is  limited  to  vehicles  of  at  most  four  different  capa- 
cities. The  total  number  of  trips  is  limited  to  50.  If  more  than  50  trips 
are  required,  data  will  be  overwritten  and  improper  and  unpredictable  function- 
ing will  result,  but  no  error  message  is  given. 


4.  OUTPUT  MAP 


From  0 through  10  output  maps  may  be  specified  on  map-bounds  cards.  If 
no  map-bounds  cards  are  present,  one  30-  by  30-inch  map  of  the  entire  region 
is  plotted.  If  more  than  10  map-bounds  cards  are  included  in  the  data  deck, 
the  cards  after  the  tenth  card  will  be  ignored.  No  message  is  printed,  and 
the  program  functioning  is  not  otherwise  affected. 


SECTION  VII 

WARNING  MESSAGE  AND  CORRECTIVE  ACTION 


One  warning  message  is  printed  by  subroutine  SECTION.  The  message  TRUCK 
CONFIGURATION  HAS  BEEN  EXTENDED  is  printed  as  the  last  line  of  the  first  page 
of  printed  output  when  more  trips  are  needed  than  either  the  minimum  number 
computed  by  PHASE2  or  the  number  specified  by  the  user.  The  problem  has  two 
possible  causes.  First,  the  time  limit  on  a trip  may  have  caused  a section 
to  be  completed  before  the  vehicle  was  filled.  In  this  case,  the  user  must 
decide  whether  the  time  limit  is  important  enough  to  necessitate  arc  extension 
of  the  number  of  trips.  Second,  the  number  of  vehicles  needed  may  have  been 
extended  because  all  segments  considered  for  addition  to  a section  contained 
enough  refuse  to  cause  the  vehicle  capacity  to  be  exceeded. 

The  problem  can  be  solved  in  one  of  two  ways.  The  simplest  requires  re- 
running PHASE2  with  a different  choice  of  starting  base  segment.  It  may  be 
necessary  to  do  this  several  times  before  the  program  runs  without  extending 
the  number  of  trips.  An  alternative  approach  is  to  cause  the  segments  that 
could  not  be  added  to  the  section  to  be  collected  one  side  at  a time.  The 
section  involved  can  be  found  by  looking  at  the  summary  of  vehicle  loads;  the 
section  (other  than  the  last)  which  has  a vehicle  load  considerably  below  the 
vehicle  capacity  will  be  the  section  in  which  the  problem  has  occurred.  Seg- 
ments near  the  last  segment  added  to  this  section  are  those  which  should  be 

modified.  The  modification  will  cause  one  or  more  of  the  segments  to  be  ser- 
viced by  two  collection  vehicles,  one  on  each  side  of  the  street. 

Note:  program  malfunction  was  observed  on  one  run.  Three  segments 

were  assigned  twice,  and  three  segments  were  not  assigned  to  any  section.  The 
cause  of  this  problem  is  still  unknown,  but  the  problem  can  be  bypassed  by 
choosing  a different  starting  base  segment. 
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SECTION  VIII 
RECOMMENDED  CHANGES 


Three  changes  in  program  PHASE2  are  recommended.  The  problem  title  could 
be  appended  to  the  section  maps.  If  the  house  count  in  array  NH  is  made 
floating-point  instead  of  integer  in  program  RCINPT,  the  array  should  be  de- 
clared type  REAL  in  the  main  program.  It  will  be  necessary  to  change  the  ab- 
solute value  function  name  from  IABS  to  ABS  in  statement  15  of  the  main  program. 
An  error  message  should  be  added  to  warn  the  user  if  more  than  50  trips  will  be 
required. 
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APPENDIX  C 


DEFINITIONS  OF  IMPORTANT  VARIABLES 


Subroutine  SHLSRT 
Subroutine  SIFTUP 
Subroutine  SORTK 
Function  IFIND 
Subroutine  NUMBER 
Subroutine  SHAPCOM 
Subroutine  COORD 
Subroutine  MAPPLT 
Subroutine  BUILD 
Subroutine  SECTION 
Program  PHASE2 


Note:  A single  variable  symbol  may  have  different  meanings  in  relation  to 

the  various  subroutines.  For  this  reason,  variables  are  defined  below 
for  each  subroutine  and  for  program  PHASE2. 


SUBROUTINE  SHLSRT 

A Array  reordered  as  array  X is  sorted 

NW  Number  of  words  to  be  sorted 

X Array  sorted  into  decreasing  order 

SUBROUTINE  SIFTUP 
L Number  of  words  in  tree 

N Pointer  to  root  at  which  ordering  of  root  with  respect  to  branches 

begins 

TREE  Array  of  distances  to  be  sorted 

SUBROUTINE  SORTK 

KN  Number  of  sorted  items  to  be  returned 

N Number  of  words  in  array  TREE 

NN  Pointer  to  root  at  which  subroutine  SIFTUP  begins  ordering  root 

with  respect  to  branches 

TREE  Array  containing  packed  distances  between  segments  and  segment 
numbers 


V 


FUNCTION  IFIND 

IARRAY  Array  being  searched 

LEN  Length  of  IARRAY 

NUM  Number  being  sought 

SUBROUTINE  NUMBER 

FORM  Output  format  for  number 

NUM  Number  to  be  plotted 

TEXT  Character  representation  of  number 
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SUBROUTINE  SHAPCOM 


AVMD 

Map  distance  conversion  factor,  in  miles  per 

map  coordinate  unit 

BR1 

Distance  to  first  break  in  segment  shape,  in 

miles 

BR2 

Distance  to  second  break  in  segment  shape,  in  miles 

ISF 

Shape  code  when  in  character  form 

R 

Radius  of  curvature  of  circular  segments,  in 

miles 

RPR 

Reciprocal  of  radius  of  curvature 

SF 

Shape  code  when  in  binary  form 

THETA 

Slope  of  line  from  starting  to  ending  node. 

in  radians 

TOTLEN 

Total  length  of  segment,  in  miles 

XNF 

X-coordinate  of  ending  node 

XNI 

X-coordinate  of  starting  node 

YNF 

Y-coordinate  of  ending  node 

YNI 

Y-coordinate  of  starting  node 

SUBROUTINE  COORD 

BR1  Distance  to  first  break  in  segment  shape,  in  miles 

BR2  Distance  to  second  break  in  segment  shape,  in  miles 

CUMLEN  Cumulative  length  along  segment,  in  miles 

RPR  Reciprocal  of  radius  of  curvature  of  a circular  segment 

S Distance  along  segment  since  previous  break 

SF  Shape  code 

XNF  X-coordinate  of  ending  node 

XNl  X-coordinate  of  starting  node 

YNF  Y-coordinate  of  ending  node 

YNI  Y-coordinate  of  starting  node 

SUBROUTINE  MAPPLT 

AVMD  Map-distance  conversion,  in  miles  per  map  coordinate  unit 

CUMLEN  Cumulative  street  length,  in  miles 

FLEN  Array  of  segment  lengths,  in  miles 

INB  Point  within  map-bounds  indicator 

ISEG  Array  of  section  assignments 


1 23 


KNODES 

NMAP 

NMAPO 


NODNUM 

NPPSEG 


SVAV 

TOTLEN 


BDATA 


ISTPR 


MINFR 


NNTEMP 

STRING 


SUBROUTINE  MAPPLT  (Concl'd.) 

Shape  code  when  in  character  form 
Count  of  nodes 

Map  strip  number  of  current  point 

Map  strip  number  of  previous  point 

Array  of  starting  node  numbers 

Array  of  ending  node  numbers 

Array  of  node  numbers 

Number  of  points  plotted  per  segment 

Height  of  map  strip,  in  inches 

Total  length  of  all  plot  strips,  in  inches 

Shape  code  when  in  binary  form 

Array  of  map  distance  conversion  factors 

Total  segment  length,  in  miles 

Array  of  node  x-coordi nates 

Array  of  node  y-coordinates 

Height  of  map  output  strips,  in  map  coordinate  units 


SUBROUTINE  BUILD 

Array  of  single  1-bit  masks 

Array  of  near-neighbor  data 

Distance  in  miles  between  street  midpoints 

Array  of  segment  numbers 

Number  of  near  neighbors  to  be  found  for  each  segment 

Number  of  words  required  to  save  near-neighbor  indicators  for  each 
segment 

Array  of  refuse  quantity,  total  traversal  time,  and  number  of  houses 
on  segments 

Number  of  segments 

Array  of  near-neighbor  segment  numbers 
Array  of  segment  numbers 

Array  of  segment  number,  refuse  quantity,  total  traversal  time,  num- 
ber of  houses,  and  near-neighbor  indicators  for  a segment 

Array  of  packed  distances  between  segments  and  segment  numbers 

Array  of  unit  (file)  numbers 
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SUBROUTINE  BUILD  (Concl'd.) 


X 

Array  of 

XT 

Array  of 

Y 

Array  of 

YT 

Array  of 

segment  midpoint  x-coordinates 
segment  midpoint  x-coordinates 
segment  midpoint  y-coordi nates 
segment  midpoint  y-coordinates 


SUBROUTINE  SECTION 

BASE  Array  of  segment  number,  refuse  quantity,  traversal  time,  number  of 
of  houses,  and  near-neighbor  indicators  for  the  current  base  segment 

BCOMP  Array  of  base  segment  near-neighbor  indicator? 

BDATA  Array  of  single  1-bit  masks 

BMINF  Array  of  base  segment  refuse  quantity,  traversal  time,  and  number  of 

houses 

COMP  Array  of  segment  near-neighbor  indicators 
CUML  Refuse  quantity  of  all  unassigned  segments 
CUMT  Traversal  time  of  all  unassigned  segments 
CURL  Refuse  quantity  on  current  segment 
CURT  Traversal  time  of  current  segment 

TRACT  Ratio  of  total  refuse  quantity  to  total  vehicle  capacity 

HISTO  Array  of  number  of  occurrences  of  number  of  shared  near  neighbors 

equal  to  HISTO  subscript 

ICODE  New  section  indicator:  ICOOE  = 0 if  a section  is  incomplete  or 
ICODE  = 1 if  a section  is  complete 

ION  Pointer  to  unassigned  segment  sharing  most  near  neighbors  with  base 

segment 

ISTO  Array  of  segment  numbers 

IST1  Array  of  counts  of  shared  near  neighbors 
IST2  Array  of  pointers  to  segment  number 
IST4  Array  of  counts  of  shared  near  neighbors 
ITR  Number  of  vehicles  of  a particular  capacity 

JON  Count  of  segments  which  share  near  neighbors  with  the  base  segment 

K Maximum  number  of  near  neighbors  found  for  any  segment 

KCUTOF  Limit  on  number  of  segments  sharing  the  most  near  neighbors  with  a 

base  segment  to  be  examined  for  inclusion  in  a section  with  the  base 
segment 

KL  Number  of  unassigned  segments 

KP  Number  of  words  required  for  near-neighbor  indicators 
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KPB 

MA 

MINFO 

N 

NBASE 

NEXTN 

NN 

NNTS 

NP 

NSTD 

NSTO 

NTO 

NTRUCK 

NUT 

PC 

PK 

SECTN 

SMLD 

STRING 

TDUMP 

TOTL 

TOTT 

TRUCK 

TRUCKS 

UNOT 

X 

Y 


SUBROUTINE  SECTION  (Concl'd.) 

Count  of  words  in  use  in  array  STRING 

Maximum  number  of  segments  in  first  section  to  be  saved  for  use  as 
base  segments 

Array  of  refuse  quantity,  traversal  time,  and  number  of  houses  on  a 
segment 

Count  of  segments 

Segment  number  of  first  base  segment  for  first  section 

Pointer  to  next  base  segment  from  NNTS  array 

Count  of  unassigned  segments 

Array  of  data  for  base  segments 

Count  of  unassigned  segments 

Count  of  base  segments  currently  saved  in  NNTS  array 

Maximum  number  of  base  segments  which  can  be  saved  in  the  NNTS  array 

Original  number  of  sections  required 

Original  number  of  vehicles  required 

Count  of  unassigned  segments  on  file  IUNIT 

Count  of  segments  assigned  to  current  section 

Count  of  segments  assigned  to  all  completed  sections 

Number  of  current  section 

Minimum  total  refuse  which  must  be  assigned  before  current  section 
is  completed 

Array  of  segment  number  followed  by  MINFO  and  COMP  arrays 
Unloading  time  at  the  landfill,  in  minutes 
Total  refuse  quantity  for  all  segments 
Total  traversal  time  for  all  segments 

Array  of  vehicle  capacity,  maximum  trip  time,  load,  actual  trip  time, 
pointer  to  first  segment  in  section,  count  of  segments  section,  and 
count  of  houses  in  section,  for  each  section 

Array  of  quantity,  capacity,  and  maximum  trip  time  for  vehicles  of 
each  capacity 

Array  of  unit  (file)  numbers 

Array  of  segment  midpoint  x-coordinates 

Array  of  segment  midpoint  y-coordinates 
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PROGRAM  PHASE  2 


FLEN 

FMPH 

1SECTN 

KCUTOF 

KN 

RNODES 

KP 

MA 

MAPS 

MINFR 

NA 

NBASE 

Ml 

NSEG 

NT 

REF 

RQF 

TC 

TMXTR 

TOTL 

TRUCK 

TRUCKS 

UNOT 

X 

XNOD 

Y 

YNOD 


Array  of  street  segment  lengths,  in  miles 
Array  of  speed  limits,  in  mph 
Count  of  sections 

Limit  on  number  of  segments  sharing  the  most  near  neighbors  with  a 
base  segment  to  be  examined  for  inclusion  in  a section  with  the  base 
segment 

Count  of  near  neighbors  to  be  found  for  each  segment 
Count  of  nodes 

Count  of  words  required  to  save  near-neighbor  indicators 

Maximum  number  of  segments  in  first  section  to  be  saved  for  use  as 
base  segments 

Count  of  section  maps  to  be  plotted 

Array  of  refuse  quantities,  total  traversal  times,  and  number  of 
houses  on  segments 

Count  of  segments 

Segment  number  of  first  base  segment  for  first  section 
Array  of  count  of  houses  on  a segment 
Count  of  segments 

Array  of  count  of  vehicles  of  each  capacity 
Refuse  quantity 

Refuse  quantity  adjustment  factor 
Array  of  vehicle  capacities 
Maximum  trip  time 
Total  refuse  quantity 

Array  of  vehicle  capacity,  maximum  trip  time,  load,  actual  trip  time, 
pointer  to  first  segment  in  section,  count  of  segments  in  section, 
and  count  of  houses  in  section,  for  each  section 

Array  of  quantity,  capacity,  and  maximum  trip  time  for  vehicles  of 
each  capacity 

Array  of  unit  (file)  numbers 

Array  of  x-coordi nates  of  segment  midpoints 

Array  of  node  x-coordinates 

Array  of  y-coordinates  of  segment  midpoints 

Array  of  node  y-coordinates 
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APPENDIX  D 


SAMPLE  PRINTED  OUTPUT 
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A 


GLOSSARY 

Air  Force  Refuse-Collection  Scheduling  Program:  a set  of  four  computer  programs 
that  perform  residential  refuse-collection  scheduling  and  produce 
printed  schedules  and  maps  of  the  routes. 

base  segment:  a segment  used  to  limit  the  addition  of  other  segments  to  its 
section  on  the  basis  of  the  number  of  neighboring  segments  common 
to  both. 

binary  search:  a procedure  for  finding  one  item  in  an  ordered  group  by  re- 
peatedly halving  the  portion  of  the  group  that  contains  the  item. 

map  coordinate  unit  (MCU):  the  length,  in  inches,  between  integral  divisions 
on  the  coordinate  system  appended  to  a map. 

node:  a numbered  point  on  a street  at  which  some  characteristic  of  the 
street  changes. 

pointer:  a variable  that  gives  the  location  of  some  other  variable. 

segment:  a portion  of  a street  between  two  nodes. 

shape  code:  characters--either  two  letters  or  a letter  followed  by  a number-- 
that  indicate  the  shape  of  a street  segment. 

spatial  clustering  of  streets:  selection  of  streets  traversed  by  a vehicle  on 
one  trip  so  that  the  streets  are  connected  by  other  streets  that 
must  be  traversed. 
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INITIAL  DISTRIBUTION 
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